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Abstract:  This  document  describes  a  GIS  approach  that  can  quickly  and 
inexpensively  identify  areas  around  military  installations  that  are 
attractive  to  residential  development  -  based  on  proposed  regional 
investments  and  policies.  Analysis  results  can  then  be  used  to  inform 
regional  planning  and  to  help  identity  long-term  impacts  of  residential 
development  on  an  installation’s  future  opportunity  to  train  and  test. 
Proposed  state  and  local  investment  in  regional  planning  policies  and 
projects  can  alter  the  expected  patterns  of  growth  and  should  be  evaluated 
with  respect  to  their  short  and  long-term  impacts  on  the  ability  to  conduct 
training  and  testing  at  the  installation.  For  demonstration  purposes,  the 
approach  is  demonstrated  for  a  multi-county  area  around  Fort  Bragg,  NC. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes.  Citation 
of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products.  All  product 
names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 

DESTROY  THIS  REPORT  WHEN  NO  LONGER  NEEDED.  DO  NOT  RETURN  IT  TO  THE  ORIGINATOR. 
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1  Introduction 

1.1  Background 

The  growth  on  the  U.S.  population  over  the  past  several  decades,  com¬ 
bined  with  advances  in  communication  and  transportation  technologies, 
has  supported  more  diffused  urban  development  patterns.  Installations 
originally  sited  in  rural  areas  are  often  now  at  the  fringe  or  in  the  midst  of 
large  urbanized  or  urbanizing  areas.  Many  of  the  increasing  restrictions  on 
training  and  testing  are  caused  indirectly  by  patterns  of  urban  develop¬ 
ment  that  result  in  incompatible  land  uses  on  either  side  of  installation 
fence  lines.  Increased  complaints  of  dust,  noise,  smoke,  and  radio  and  TV 
interference  are  directly  related  to  the  development  of  urban  patterns  near 
installations.  Increases  in  pressure  to  protect  endangered  species  and  sen¬ 
sitive  habitat  is  directly  related  to  the  percent  of  the  habitat  in  the  region 
on  the  military  installations.  Appropriately  dark  training  areas  can  be  lost 
as  city  lights  follow  urban  development. 

A  tool  that  could  predict  patterns  of  urban  development  would  be  valuable 
to  help  installations  better  plan  for,  or  entirely  avoid  such  land  use  con¬ 
flicts.  Various  techniques  have  been  developed  and  applied  to  predict  the 
potential  for  the  incompatible  land  uses  around  military  installations, 
ranging  from  simple  questionnaires  answered  by  installation  personnel  to 
more  complex,  spatially  explicit  dynamic  simulation  modeling. 

Perhaps  the  most  common  approach  involves  identifying  1-,  3-,  and/or  5- 
mile  buffers  around  installations  using  current  and  historic  digital  land 
use  maps,  and  then  by  measuring  the  urbanized  area  in  these  areas  using  a 
time-series  of  maps.  This  method  creates  a  simple  graph  of  time  vs.  total 
amount  of  urbanized  area  (Lozar  2003).  These  simple  trends  can  be  very 
useful  and  often  sufficient  to  predict  growth  based  on  past  trends.  How¬ 
ever,  extrapolating  future  growth  based  on  past  trends  can  be  misleading 
(Figure  1).  Growth  might  continue  to  accelerate  where  plenty  of  land  is 
available  for  development  and  where  the  fringe  of  a  major  urbanized  area 
represents  an  area  of  economic  interest.  Growth  might  moderate  if  an  ur¬ 
banized  center  is  already  in  the  buffer  area,  and  is  growing  along  the  edge 
of  an  installation.  Growth  might  decelerate  if  the  there  is  little  or  no  more 
land  available  to  develop. 
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Figure  1.  Past  trends  are  not  necessarily  good  at  predicting  the  future. 

A  more  accurate  prediction  of  urban  development  patterns  (and  potentials 
for  incompatible  land  use)  would  require  a  careful  analysis  of  the  spatial 
relationships  of  growing  urban  centers,  available  land,  natural  barriers 
(e.g.,  rivers)  or  man-made  barriers  (e.g.,  limited  access  highways),  and 
zoning.  To  date,  the  more  involved  procedures  are  very  expensive,  time 
consuming,  and  have  substantial  data  requirements.  Consequently,  there 
is  a  need  to  describe  and  demonstrate  a  quick  and  inexpensive  GIS-based 
analysis  that  uses  nationally  available  digital  maps  to  identify  the  relative 
attractiveness  of  land  in  the  vicinity  of  an  installation  to  residential  devel¬ 
opment  in  response  to  proposed  regional  planning  policies  and  invest¬ 
ments.  The  results  would  then  be  useful  to  explain  the  long-term  direct 
and  indirect  consequences  of  proposed  regional  investments  and  policies 
on  local  military  training  and  testing  opportunities. 

1.2  Objectives 

The  objective  of  this  study  was  to  create  an  inexpensive  GIS-based  ap¬ 
proach  for  predicting  future  urban  growth  patterns  that: 

1.  Uses  nationally  available  data 

2.  Requires  minimal  data  preparation 

3.  Is  readily  calibrated 

4.  Does  not  require  time-series  data 

5.  Identifies  the  relative  attractiveness  of  areas  around  military  installa¬ 
tions  much  better  than  the  current  image-processing  based  trend 
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analyses  and  much  less  expensively  than  currently  available  urban 
growth  simulation  models. 

1.3  Approach 

The  analysis  is  first  described  in  detail  and  then  a  sample  application  is 
explored.  All  instructions  and  scripts  are  included  that  will  allow  a  GIS 
specialist  familiar  with  UNIX  to  conduct  the  analysis  for  their  areas  of  in¬ 
terest. 

1.4  Scope 

Note  that  this  work  describes  an  approach  for  predicting  future  urban 
growth  patterns.  Results  generated  by  the  approach  are  not  ready  to  be 
used  in  informing  local  planning  processes. 

Many  of  the  encroachment  factors  affecting  training  are  a  direct  result  of 
the  regional  urban  pattern.  This  pattern  is  the  result  in  many  areas  of  free- 
market  property  exchange  and  ownership.  Land  attractiveness  to  the  free 
market  participants  (e.g.,  homeowners)  is  based  on  local  county  and  city 
investments  and  policies  including  roads,  utilities,  and  zoning.  Future 
training  and  testing  capacities  on  installations  will  be  directly  and  indi¬ 
rectly  affected  by  regional  policies  and  investments  so  that  those  capacities 
are  also,  in  part,  a  function  of  regional  land  use  investments  and  policies 
(Figure  2). 
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Figure  2.  Cause-effect  chains. 
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The  effort  described  in  this  report  produced  a  map  showing  the  relative 
attractiveness  of  properties  across  the  landscape  to  future  urbanization. 
The  model  created  to  generate  this  map  can  then  be  used  to  test  for  alter¬ 
native  policies  and  investments  including  zoning,  property  purchases,  and 
investments  in  roads. 

1.5  Mode  of  Technology  Transfer 

It  is  anticipated  that  the  information  in  this  report  and  accompanying  digi¬ 
tal  maps  will  directly  assist  military  installations,  their  supporting  organi¬ 
zations,  and  the  Office  of  Economic  Adjustment  (OEA).  This  report  will  be 
made  accessible  through  the  World  Wide  Web  (WWW)  at  URL: 

http :  //www.  cecer.  a  rm  y.  m  i  I 
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2  Review  of  Urban  Growth  Models 

Different  organizations  have  developed  many  urban  growth  simulation 
models  for  a  variety  of  purposes  ranging  from  research  and  education  to 
city  and  regional  planning  (USEPA  2000).  Each  model  captures  a  different 
combination  of  the  factors  involved  in  growth  to  answer  specific  questions. 
Generally,  models  can  be  separated  into  two  distinct  spatial  scales:  city 
and  region.  Models  such  as  DRAM/EMPAL  and  MEPLAN  are  popular  in 
both  the  United  States  and  overseas,  and  focus  on  identifying  growth  by 
income  and  housing  costs.  These  and  other  models  focus  on  the  city  itself 
and  deal  with  growth  over  the  short  term  (no  more  than  a  couple  of  dec¬ 
ades). 

Regional  models  address  questions  regarding  gross  growth  patterns  over 
the  long  term  (many  decades)  and  the  associated  implications  with  respect 
to  environmental  measures.  Regional  models  can  help  military  installa¬ 
tions  address  such  questions  as: 

•  Where  are  people  likely  to  be  living  in  the  next  50  years? 

•  What  percent  of  important  habitat  in  the  region  will  be  on  local  mili¬ 
tary  installations? 

•  What  limitations  will  regional  land  use  patterns  place  on  military  in¬ 
stallation  activities? 

Models  that  specifically  address  regional  growth  patterns  over  the  course 
of  many  decades  include  the  California  Urban  Futures  model,  version  2 
(CUF-2),  SLEUTH,  Landuse  Evolution  Assessment  Model  (LEAM™), 
Smart  Places,  and  What  If?: 

•  CUF-2  uses  a  set  of  econometric  models  to  project  future  population, 
household,  and  employment.  The  landscape  is  gridded  into  one- 
hectare  developable  land  units  (DLU).  The  value  of  each  DLU  for  each 
competing  land  use  is  computed  and  land  use  change  is  based  on  com¬ 
paring  these  relative  prices. 

•  SLEUTH  derives  its  name  from  the  types  of  raster  GIS  data  inputs: 
slope,  land  use,  urban,  exclusion,  transportation,  and  hillshading 
(Clarke  and  Gaydos  1998).  SLEUTH  is  a  cellular  automaton  model  that 
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updates  the  state  of  each  cell  in  each  time  step  based  on  the  state  of 
each  cell  and  its  immediate  neighbors. 

•  LEAM™  takes  a  similar  approach  (Deal  2003).  LEAM  uses  a  30-meter 
grid-cell  cellular  automation  approach,  in  which  each  cell  is  assigned  a 
variety  of  development  attractor  values  and  the  landscape  evolves 
based  on  equations  that  assign  future  cell  states  based  on  the  state  of 
cells  and  their  immediate  neighbors. 

•  Smart  Places  is  implemented  as  an  extension  to  ESRI’s  ArcView  GIS, 
and  was  designed  to  be  easy  to  use  for  those  familiar  with  that  system 
(ESRI 1999). 

•  What  If?  is  a  stand-alone  product  that  uses  ESRI  Shapefiles  (Kloster- 
man  2003).  Parcels  are  represented  as  entities  using  the  ArcView 
Shapefiles  that  change  through  time  based  on  the  growth  needs  of  the 
region,  the  state  of  the  parcel,  and  the  state  of  its  neighbors. 

To  be  useful  for  projection,  all  urban  growth  models  must  be  calibrated  (a 
process  that  can  be  time  consuming  and  difficult).  The  SLEUTH  model  has 
been  calibrated  by  testing  thousands  of  combinations  of  coefficients  using 
very  fast  computers.  LEAM  model  was  calibrated  using  a  statistical  tech¬ 
nique  (George  Gertner,  personal  communication).  These  and  other  models 
can  be  run  for  a  location  using  default  calibrations,  which  can  be  useful  for 
educational  purposes  and  for  generating  thoughtful  conversation  among 
stakeholders  and  other  participants. 

Many  land  use  change  models  have  been  developed  and  are  being  used  to 
test  alternative  land  use  policies  with  respect  to  their  impact  on  future 
land  patterns  in  and  around  cities  and  towns  (USEPA  2000).  The  Corps  of 
Engineers  is  adopting  LEAM™  to  help  evaluate  how  alternative  regional 
policies  and  land  ownership  patterns  affect  future  land  development.  The 
Corps’  primary  interest  is  to  help  minimize  future  land  use  conflict  result¬ 
ing  from  the  development  of  new  uses  in  areas  that  are  and  will  be  im¬ 
pacted  by  military  training  and  testing  activities. 

Although  prediction  of  urban  growth  is  more  art  than  science,  it  is  ex¬ 
tremely  important  as  part  of  the  necessary  conversations  among  stake¬ 
holders.  In  the  short  term  (e.g.,  up  to  10-15  years),  urban  development  will 
likely  proceed  along  the  lines  established  by  stakeholders— especially 
landowners  and  those  who  invest  in  anticipation  of  development.  Such 
stakeholders  follow  strategic  purchases  of  land  with  aggressive  and  persis- 
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tent  pressure  on  and  participation  in  the  political  processes  concerned 
with  infrastructure  investments  and  zoning  issues.  Urban  growth  models 
developed  to  predict  the  short-term  patterns  must  capture  the  goals,  ambi¬ 
tions,  and  expectations  of  landowners,  investors,  and  the  political  proc¬ 
esses.  Modeling  at  this  scale  will  be  extremely  data  intensive;  and  the 
model  itself  can  actually  change  the  process  and  therefore  the  future. 

In  the  long  term  (30,  50,  and  more  years),  unanticipated  changes  in  tech¬ 
nology  can  dramatically  change  the  way  land  use  patterns  develop.  For  ex¬ 
ample,  the  advent  of  the  automobile  allowed  unprecedented  connections 
of  cities  to  surrounding  agricultural  land.  Adding  tractors  and  other  auto¬ 
mated  implements  to  the  farm  continues  to  reduce  the  farm  labor  re¬ 
quirement.  The  percent  of  population  in  the  cities  has  consequently 
jumped  from  6  percent  in  1800  to  40  percent  in  1900  and  75  percent  in 
2000  (Figure  3). 

Future  technology  drivers  may  continue  this  trend  or  begin  reversing  it,  or 
change  the  patterns  of  human  settlement  in  entirely  new  ways.  For  exam¬ 
ple,  businesses  that  rely  primarily  on  information  exchange  may  soon  have 
sufficient  Internet  bandwidth  to  support  remote  employees  to  a  far  greater 
extent.  A  home  office  may  integrate  a  wall  “window”  that  can  be  virtually 
shared  with  anyone,  anywhere  with  a  similar  device.  With  a  simple  flip  of  a 
switch,  people  may  virtually  share  office  space  across  the  world,  as  needed, 
nearly  eliminating  the  need  for  businesses  to  support  a  central  office  loca¬ 
tion.  When  the  “window”  is  not  being  used  for  business,  the  office  worker 
may  switch  to  any  of  hundreds  of  live  scenes  from  across  the  world,  the 
moon,  or  Mars.  With  no  need  to  physically  “go  to  work,”  the  rural  land¬ 
scape  may  see  a  dramatic  infusion  of  business  people. 
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Time 


Figure  3.  Rural  to  urban  population  shift. 

Urban  growth  projection  models  must  be  applied  very  carefully.  The  po¬ 
tential  for  a  single  model  to  predict  urban  growth  anywhere  is  very  limited. 
Differences  in  physical  and  economic  geography  make  it  necessary  to  cali¬ 
brate  any  model  for  the  local  circumstances.  Differences  in  culture,  demo¬ 
graphic  patterns,  and  personal  income  levels  will  result  in  different  overall 
settlement  patterns.  For  example,  the  poor  tend  to  live  more  densely  and 
simply,  while  the  rich  tend  to  own  and  occupy  land  at  much  lower  densi¬ 
ties.  Calibration  must  account  for  current,  historic,  and  anticipated  demo¬ 
graphic  makeup  of  populations. 
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3  Approach  for  Predicting  Future  Urban 
Patterns 

3.1  Overview 

Urban  development  in  the  United  States  results  from  interactions  among 
political  forces,  commercial  development,  and  individual  homeowners.  In¬ 
vestments  through  political  processes  into  various  infrastructure  projects 
including,  roads,  highways,  interstate  highways,  and  bridges.  The  same 
processes  develop  policies  in  the  form  of  zoning,  tax  incentives,  and  ex¬ 
change  of  property  rights.  Commercial  developers  create  neighborhoods, 
industrial  parks,  and  shopping  areas.  Individual  homebuyers  purchase 
homes.  People  make  all  the  associated  decisions,  which  they  base  on  for¬ 
mal  and  informal  analyses  of  the  attractiveness  of  locations  for  meeting 
various  goals  and  objectives.  The  challenge  of  this  study  is  to  quickly  and 
inexpensively  convert  readily  available  GIS  maps  into  projections  of  where 
people  will  live  (Figure  4).  The  following  steps  identity  the  attractiveness 
of  undeveloped  land  to  residential  development: 

1.  Acquire  land  use  data  for  the  area  of  interest. 

2.  Identify  a  list  of  attractors  for  growth  (human  habitat  suitability). 

3.  For  each  attractor: 

a.  develop  the  attractor  map 

b.  develop  a  probability  of  finding  residential  areas  within  each  attrac¬ 
tor  category 

c.  convert  the  attractor  map  to  a  probability  map. 

4.  Combine  all  probability  maps  into  an  overall  attractor  map. 

5.  Develop  a  probability  of  finding  residential  areas  within  each  attractor 
category. 

6.  Convert  the  attractor  map  to  a  probability  map. 

The  prediction  of  future  land  development  is  based  on  the  demonstration 
of  land  preferences  for  residential  development  based  on  historic  decisions 
reflected  in  the  current  land  use  patterns.  The  U.S.  Geological  Survey 
(USGS)  National  Land  Cover  Data  (NLCD)  program  provides  a  standard 
set  of  readily  available  land  use  data  across  the  entire  United  States  in  the 
early  1990s  and  a  partial  set  representing  the  early  2000s.  Data  from  ei¬ 
ther  time  period  is  required. 
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Figure  4.  The  maps  on  the  left  are  to  be  turned  into  a  future  land  use  map. 


The  next  step  is  to  develop  a  list  of  residential  attractors  that  are  simple  to 
address,  sufficient  to  predict  urban  patterns,  and  sufficient  across  socio¬ 
economic  classes.  (Note  that  the  NLCD  data  is  imagery  based  and  makes 
no  distinction  between  socio-economic  classes.)  For  the  demonstrations 
reported  here,  the  factors  considered  are: 

•  current  land  use 

•  slope 

•  distances  to: 

a.  water  and  forest 

b.  cities  (city  centers  /  employment  centers) 

c.  county  road 

d.  state/Federal  highway 

e.  limited  access  highway 

f.  nearest  intersection. 


These  factors  are  considered  without  specific  regard  for  socioeconomic 
level  because,  to  predict  potential  change  in  training/testing  capacity  at  an 
installation  and  in  its  surrounding  region,  it  is  most  important  to  predict 
locations  of  residences,  which  can  be  sensitive  to  military  training  activity. 
Also,  while  there  is  ready  availability  of  input  maps  that  indicate  location 
of  residential  areas,  acquiring  data  indicating  the  socioeconomic  level  of 
neighborhoods  can  be  expensive. 

The  factors  listed  above  were  chosen  for  their:  (1)  apparent  importance  in 
consideration  of  constructing  residences,  and  (2)  universal  availability  of 
supporting  information  across  the  United  States.  Factors  other  than  those 
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in  the  above  list  can  be  considered  by  acquiring  or  developing  additional 
input  maps  using  the  algorithm  described  below.  Next,  the  attractor  maps 
must  be  developed  in  three  general  steps: 

1.  Acquisition  of  base  data 

2.  Manual  processing  of  GIS  data 

3.  Machine  processing  of  data. 

Ideally,  acquisition  relies  on  no  more  than  downloading  nationally  avail¬ 
able  maps  from  public  (or  at  least  inexpensive)  web  sites.  Manual  process¬ 
ing  must  be  minimized  to  allow  quick  application  to  multiple  sites.  Typi¬ 
cally,  manual  steps  that  are  difficult  to  automate  include  assembling  the 
base  data,  reformatting  the  data  into  a  common  format,  and  converting 
the  coordinate  system  and  projections  to  a  selected  standard. 

Finally,  processes  that  can  be  automated  to  create  the  required  attractor 
maps  are  applied  to  the  results  of  the  manual  processing.  The  attractor 
maps  are  factually  based  and  represent  information  in  relevant  units,  such 
as  miles,  meters,  slope  percent,  and  minutes.  For  example,  an  attraction- 
to-cities  map  will  be  expressed  in  minutes  of  travel  time;  an  attraction-to- 
water  map  may  be  shown  in  meters. 

Finally,  each  attractor  map  is  individually  converted  to  a  residential- 
probability  map.  This  is  done  by  developing  a  graph  relating  the  percent  of 
land  development  to  the  attractor  map  categories  (Figure  5).  In  the  graph, 
the  Data  Flows  from  left  to  right  following  the  arrows,  calculating  the 
probability  of  residential  development  from  the  input  maps  on  the  left. 

The  input  maps  are  used  in  different  combinations  to  create  a  number  of 
interim  maps,  which  are  used  together  to  calculate  the  residential  prob¬ 
ability  map  (Figure  6).  This  graph  is  developed  by  cross-tabulating  the 
number  of  residential  and  potential  residential  locations  with  attractor 
ranges.  The  full  range  of  the  attractor  map  values  is  divided  into  sub¬ 
ranges  (e.g.,  10).  The  percent  developed  in  each  range  is  calculated  as  the 
number  of  residential  cells  divided  by  the  sum  of  the  residential  and  po¬ 
tentially  residential  cells.  Residential  cells  are  NLCD  categories  21  and  22 
and  potential  residential  cells  are  all  other  categories  that  could  be  devel¬ 
oped  into  residential  cells  (avoiding  categories  like  water,  swamp,  high¬ 
ways,  etc.) 


P' 
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Figure  6.  Graph  relating  an  attractor  to  probability  of  residential 

development. 

A  wide  variety  of  graphs  are  possible.  For  example,  the  top  left  graph  in 
Figure  7  might  represent  the  probability  of  development  surrounding 
some  negative  influence  -  like  an  oil  refinery.  The  top  right  is  typical  of 
positive  attractors  such  as  city  centers.  A  variant  shown  in  the  bottom 
right  shows  a  different  power  in  the  attractiveness  decay.  The  bottom  left 
indicates  an  attractor  that  has  little  influence  on  the  attractiveness.  Forest 
or  water  in  some  parts  of  the  country  might  have  a  strong  positive  influ¬ 
ence,  but  in  others  where  water  and  trees  are  extremely  common,  the  in¬ 
fluence  might  be  represented  by  the  lower-left  graph. 


Figure  7.  Range  of  potential  attractor  to  percent  developed  graphs. 
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At  this  point,  each  attractor  is  associated  with  a  respective  percent- 
developed  graph  (values  o-i.o)  and  each  attractor  map  can  be  converted, 
with  this  graph,  to  a  probability  map  for  that  attractor.  The  probability 
maps  are  then  integrated  into  a  combined  attractor  map  by  adding  all  the 
probability  maps,  then  dividing  the  result  by  the  total  number  of  maps 
This  results  in  an  index  map  with  values  between  o  and  l.o.  This  combined 
attractor  map  is  then  processed  in  a  manner  identically  to  the  individual 
attractor  maps:  (l)  the  combined  map  is  divided  into  a  number  of  ranges; 
(2)  the  percentage  of  land  that  could  be  residential  that  actually  is  resi¬ 
dential  is  determined;  and  (3)  a  graph  is  generated  from  that  data;  and 
(4)  it  is  used  to  convert  the  combined  attractor  map  into  a  combined  prob¬ 
ability  map.  This  is  the  final  product,  which  indicates  the  relative  attrac¬ 
tiveness  of  land  to  residential  development.  Chapter  4  details  this  process. 

Figure  8  details  the  steps  outlined  in  Figure  4;  it  shows  the  steps,  proc¬ 
esses,  and  results  of  creating  the  maps  described  in  Figure  5.  The  grey  ar¬ 
rows  represent  three  key  steps  once  starting  maps  have  been  downloaded 
from  the  Internet: 

1.  Projection  and  coordinate  systems  conversion 

2.  GIS  hand  processing 

3.  GIS  scripts. 
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Figure  8.  Computer  and  human  steps. 
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3.2  GIS  Hand  Processing 

The  description  of  each  of  the  “hand  processed”  maps  in  the  top  center  of 
Figure  8  are  described  below.  Each  must  contain  information  with  specific 
category  information  as  indicated  below.  (Information  in  Appendix  B  will 
help  GIS  technicians  develop  these  maps.) 

3.2.1  smallCity 

This  is  a  user-developed  raster  map  identifying  the  central  areas  of  small 
cities  to  represent  an  attractor  point  to  shopping  and  employment  loca¬ 
tions.  The  size  of  these  cities  is  user  defined,  but  an  average  population 
size  of  the  cities  must  also  be  determined  to  assist  in  the  attractiveness 
processing  to  these  cities.  The  Census  Bureau  maps  can  be  used  as  a  guide, 
but  it  is  the  city  business  centers  that  need  to  be  developed. 

Categories: 

o  -  no  city 
l  -  city  center  area. 

3.2.2  mediumCity 

This  is  identical  to  the  smallCity  map  except  for  medium  cities. 

Categories: 

o  -  no  city 
l  -  city  center  area. 

3.2.3  largeCity 

This  is  identical  to  the  mediumCity  map  except  for  large  cities. 

Categories: 

o  -  no  city 
l  -  city  center  area. 

3.2.4  xlCity 

This  is  identical  to  the  largeCity  map  except  for  the  largest  cities. 
Categories: 

o  -  no  city 
l  -  city  center  area. 
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3.2.5  noGrowth_att 

This  user-developed  map  identifies  where  residential  development  cannot 
occur  due  to  property  owner  preferences  or  to  local  zoning  requirements. 
No  growth  areas  are  derived  from  property  ownership  maps  including 
Federal  land,  and  from  landcover  maps. 

Categories: 

o  -  growth  permitted 
l  -  residential  development  not  allowed. 

3.2.6  dem 

This  is  a  digital  elevation  model  containing  integer  values  of  meters  above 
sea  level.  Generally,  digital  elevation  models  are  useful  for  this  application 
without  any  special  processing. 

3.2.7  boundary 

This  map  identifies  the  extent  within  which  the  map  processing  is  accom¬ 
plished  -  the  study  area.  The  extent  is  typically  defined  by  a  set  of  contigu¬ 
ous  counties. 

Categories: 

o  -  outside  the  study  area 
l  -  inside  the  study  area. 

3.2.8  landcover 

This  map  is  based  on:  “boundary”  (above) 

This  user-provided  map  can  be  based  on  the  USGS  NLCD  series  and,  in 
fact,  uses  the  NLCD  categories.  Advanced  users  might  develop  their  own 
maps  from  processed  imagery.  The  map  is  clipped  using  the  study-area 
extent  map  -  boundary. 

Categories  (NLCD): 

o  -  No  classification 

n  -  Open  water 

12  -  Perennial  Ice/Snow 

21  -  Low  Intensity  Residential 

22  -  High  Intensity  Residential 

23  -  Commercial/Industrial/Transportation 
31  -  Bare  Rock/Sand/Clay 
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32  -  Quarries/Strip  Mines/Gravel  Pits 

33  -  Transitional 

41  -  Deciduous  Forest 

42  -  Evergreen  Forest 

43  -  Mixed  Forest 
51  -  Shrubland 

61  -  Orchards /Vineyards/Other 
71  -  Grasslands/Herbaceous 

81  -  Pastures/ Hay 

82  -  Row  Crops 

83  -  Small  Grains 

84  -  Fallow 

85  -  Urban/Recreational  Grasses 

91  -  Woody  Wetlands 

92  -  Emergent  Herbaceous  Wetlands. 

3.2.9  interstates 

This  map  is  based  on:  “boundary”  (page  16) 

This  is  also  user  supplied,  but  based  on  readily  available  maps  provided 
through  the  internet.  It  is,  like  all  of  the  other  maps,  raster  based. 

Category: 

Limited  access  highway. 

3.2.10  otherroads 

This  map  is  based  on:  “boundary”  (page  16) 

This  is  also  user  supplied,  but  based  on  readily  available  maps  provided 
through  the  internet.  It  is,  like  all  of  the  other  maps,  raster  based. 

Categories: 
o  -  no  road 

1  -  Limited  access  highway 

2  -  Federal  highway 

3  -  State  highway 

4  -  County  road 
5- 

6  -  Ramps  to  limited  access  highways. 

3.3  Development  of  Attractor  Maps 

The  maps  in  this  section  and  those  following  are  generated  with  a  GIS 
script.  This  section  first  lists  the  attractor  maps  that  hold  factual  informa- 
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tion  about  the  relationship  of  each  cell  to  surrounding  attractors  (they  end 
in  _att).  The  remainder  of  the  maps  are  interim  ones,  upon  which  the  at¬ 
tractor  maps  are  based. 

3.3.1  cities_att 

This  map  is  based  on: 

smallCityTime  (page  21) 
mediumCityTime  (page  21) 
largeCityTime  (page  21) 
xlCityTime  (page  21). 

This  map  combines  the  travel  times  to  the  nearest  cities  from  the  four 
categories  into  an  overall  attraction  to  cities  map.  This  map  represents  the 
attractiveness  of  cities  based  on  their  employment,  shopping,  and  human 
contact  possibilities. 

The  steps  of  the  algorithm  consist  of: 

1.  The  average  population  of  each  city  category  is  divided  by  the  re¬ 
spective  driving  time  and  the  four  results  are  summed. 

2.  This  result  is  histogram  stretched  to  a  scale  of  0-255. 

3.3.2  slope_att 

This  map  is  based  on:  “dem”  (page  16). 

The  steps  of  the  algorithm  consist  of: 

1.  A  standard  raster-GIS  analysis  technique  is  used  that  calculates  the 
maximum  slope  from  the  center  of  each  cell  to  one  of  the  eight  sur¬ 
rounding  cells. 

2.  The  results  are  given  in  degrees. 

3.3.3  ramp_att 

This  map  is  based  on: 

travelTime30  (page  22) 
travelTimeqo  (page  23) 
interstates  (page  17) 
otherroads  (page  17). 

The  steps  of  the  algorithm  consist  of: 

1.  Find  all  roads  categorized  as  ramps  in  the  road  map  that  border 
limited  access  highways. 
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2.  Compute  the  cumulative  driving  time  from  these  out  to  a  total  driv¬ 
ing  time  of  30  minutes  at  a  30-meter  resolution  -  saving  the  re¬ 
sults. 

3.  Compute  the  cumulative  driving  time  from  the  30-minute  edge  in 
the  previous  step  to  the  remainder  of  the  study  area  at  a  90-meter 
resolution. 

4.  Combine  the  results  at  a  30-meter  resolution  with  the  first  results 
taking  precedence  over  the  second. 

3.3.4  staterd_att 

This  map  is  based  on: 

travelTime30  (page  22) 
travelTimepo  (page  23) 
interstates  (page  17) 
other  roads  (page  17). 

The  steps  of  the  algorithm  consist  of: 

1.  Find  all  roads  categorized  as  state  highways  in  the  road  map. 

2.  Compute  the  cumulative  driving  time  from  these  out  to  a  total  driv¬ 
ing  time  of  30  minutes  at  a  30-meter  resolution  -  saving  the  re¬ 
sults. 

3.  Compute  the  cumulative  driving  time  from  the  30-minute  edge  in 
the  previous  step  to  the  remainder  of  the  study  area  at  a  90-meter 
resolution. 

4.  Combine  the  results  at  a  30-meter  resolution  with  the  first  results 
taking  precedence  over  the  second. 

3.3.5  intersect_att 

This  map  is  based  on: 

travelTime30  (page  22) 
travelTimepo  (page  23) 
intersection  (page  25). 

The  steps  of  the  algorithm  consist  of: 

1.  Compute  the  cumulative  driving  time  from  intersections  out  to  a  to¬ 
tal  driving  time  of  30  minutes  at  a  30-meter  resolution  -  saving  the 
results. 

2.  Compute  the  cumulative  driving  time  from  the  30-minute  edge  in 
the  previous  step  to  the  remainder  of  the  study  area  at  a  90-meter 
resolution. 
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3.  Combine  the  results  at  a  30-meter  resolution  with  the  first  results 
taking  precedence  over  the  second. 

3.3.6  water_att 

This  map  is  based  on: 
landcover  (page  16). 

The  steps  of  the  algorithm  consist  of: 

1.  Find  all  cells  in  the  landcover  map  categorized  as  open  water  or 
swamp. 

2.  Generate  1,  2,  3,  and  4  cell  buffer  zones  around  these  cells. 

3.  Multiply  these  to  make  sure  they  represent  meter-distances  from 
cells. 

3.3.7  forest_att 

This  map  is  based  on: 
landcover  (page  16). 

The  steps  of  the  algorithm  consist  of: 

1.  Find  all  cells  in  the  landcover  map  categorized  as  forest. 

2.  Generate  1,  2,  3,  and  4  cell  buffer  zones  around  these  cells. 

3.  Multiply  these  to  make  sure  they  represent  meter-distances  from 
cells. 

3.3.8  road_att 

This  map  is  based  on: 

travelTime30  (page  22) 
travelTimepo  (page  23) 
interstates  (page  17) 
otherroads  (page  17). 

The  steps  of  the  algorithm  consist  of: 

1.  Identify  roads  classified  as  county  roads  in  the  road  map. 

2.  Compute  the  cumulative  driving  time  from  these  out  to  a  total  driv¬ 
ing  time  of  30  minutes  at  a  30-meter  resolution  -  saving  the  re¬ 
sults. 

3.  Compute  the  cumulative  driving  time  from  the  30-minute  edge  in 
the  previous  step  to  the  remainder  of  the  study  area  at  a  90-meter 
resolution. 
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4.  Combine  the  results  at  a  30-meter  resolution  with  the  first  results 
taking  precedence  over  the  second. 

3.3.9  smallCityTime 

This  map  is  based  on: 

smallCity  (page  14) 
travelTime30  (page  22) 
travelTimepo  (page  23) 

The  steps  of  the  algorithm  consist  of: 

1.  Compute  the  cumulative  driving  time  from  small  cities  out  to  a  total 
driving  time  of  30  minutes  at  a  30-meter  resolution  -  saving  the  re¬ 
sults. 

2.  Compute  the  cumulative  driving  time  from  the  30-minute  edge  in 
the  previous  step  to  the  remainder  of  the  study  area  at  a  90-meter 
resolution. 

3.  Combine  the  results  at  a  30-meter  resolution  with  the  first  results 
taking  precedence  over  the  second. 

Result: 

The  result  is  a  driving  time  map  for  every  cell  in  the  study  area  to  the 
nearest  small  city  -  with  greater  fidelity  within  for  locations  within  30- 
minute  driving  times. 

3.3.10  mediumCityTime 

This  map  is  based  on: 

mediumCity  (page  15) 
travelTime30  (page  22) 
travelTimepo  (page  23) 


The  steps  of  the  algorithm  consist  of: 

Identical  to  the  smallCityTime  algorithm,  except  with  driving  times 
computed  to  the  nearest  medium  city. 

3.3.11  largeCityTime 

This  map  is  based  on: 

largeCity  (page  15) 
travelTime30  (page  22) 
travelTimepo  (page  23) 
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The  steps  of  the  algorithm  consist  of: 

Identical  to  the  smallCityTime  algorithm,  except  with  driving  times 
computed  to  the  nearest  large  city. 

3.3.12  xlCityTime 

This  map  is  based  on: 

xlCity  (page  15) 
travelTime30  (page  22) 
travelTimepo  (page  23) 


The  steps  of  the  algorithm  consist  of: 

Identical  to  the  smallCityTime  algorithm,  except  with  driving  times 
computed  to  the  nearest  extra-large  city. 

3.3.13  travelTime90 

This  map  is  based  on:  “travelSpeedpo”  (page  23) 

The  steps  of  the  algorithm  consist  of: 

Multiply  the  travel  speed  to  get  the  equivalent  travel  time  for  traversing 
30  meters. 

3.3.14  travelTime30 

This  map  is  based  on:  “travelSpeed3o”  (page  22) 

The  steps  of  the  algorithm  consist  of: 

Multiply  the  travel  speed  to  get  the  equivalent  travel  time  for  traversing 
30  meters. 

3.3.15  travelSpeed30 

This  map  is  based  on: 

noGrowth_att  (page  16) 
boundary  (page  16) 
hwyBuff  (page  25) 
roadTravelSpeed30  (page  23) 
landTravelSpeed30  (page  24). 

The  steps  of  the  algorithm  consist  of: 

1.  Combine  the  input  maps  using  the  following  logic: 
if  the  cell  is  outside  the  study  area,  set  it  to  null 
otherwise,  if  the  cell  travel  speed  is  null,  set  it  to  null 
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otherwise  if  the  road  travel  speed  is  greater  than  o  set  it  to  the  road 
travel  speed, 

otherwise,  if  growth  is  not  allowed  in  the  cell,  set  it  to  null, 
otherwise,  if  the  cell  is  within  the  limited  access  highway  set  it  to 
null, 

otherwise,  set  it  to  the  land  travel  speed. 

2.  The  result  is  a  map  that  provides  the  travel  speed  across  every  grid 
cell  for  which  travel  is  permitted  or  possible. 

3.3.16  travelSpeed90 

This  map  is  based  on: 

hwyBuff  (page  25) 
roadTravelSpeed30  (page  23) 
landTravelSpeed30  (page  24) 
travelableLandSpeed  (page  24)1 

The  steps  of  the  algorithm  consist  of: 

1.  Although  this  results  in  values  useful  for  later  analyses  at  90-meter 
resolutions,  it  is  calculated  using  30-meter  resolution  cells. 

2.  Find  the  maximum  road  travel  speed  in  the  surrounding  9-cell 
neighborhood. 

3.  Find  the  maximum  land  travel  speed  in  the  surrounding  9-cell 
neighborhood. 

4.  Find  out  if  there  is  a  limited-access-highway  buffer  in  the  9-cell 
neighborhood. 

5.  If  the  maximum  road  speed  is  non-zero,  choose  that  speed 
otherwise,  if  there  is  a  highway  buffer,  set  to  null 
otherwise,  choose  the  maximum  overland  speed. 

3.3.17  roadTravelSpeed30 

This  map  is  based  on: 

interstates  (page  17) 
otherroads  (page  17). 

This  is  simply  an  assignment  of  travel  speeds  (in  miles  per  hour  [mph]). 
Note  that  the  travel  speed  on  roads  is  set  and  assumes  that  road  speed 
does  not  change  with  volume.  Over  the  long  term,  this  might  be  a  reason¬ 
able  assumption  if  cities  respond  to  congestion  by  widening  roads.  Travel 
speed  assignments  are  made  as  follows: 
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Cat 

Description 

Speed  (mph) 

0 

No  road 

0 

1 

Limited  access  highway 

70 

2 

Federal  highway 

50 

3 

State  highway 

40 

4 

County  road 

30 

5 

20 

6 

Ramps  to  limited  access  highways 

30 

7 

10 

3.3.18  landTravelSpeed30 

This  map  is  based  on: 

landcover  (page  16) 
noGrowth_att  (page  16). 

This  also  is  a  simple  assignment  of  overland  travel  speeds.  This  accommo¬ 
dates  the  notion  that  driveways  to  homes  can  be  developed  off-road.  No¬ 
growth  areas  are  assigned  the  null  value  to  indicate  the  inability  to  travel 
there.  Otherwise  travel  speed  assignments  are  made  based  on  the  land- 
cover  categories: 


Cat 

Description 

Speed  (mph) 

0 

No  classification 

0 

11 

Open  water 

0 

12 

Perennial  Ice/Snow 

0 

21 

Low  Intensity  Residential 

0.5 

22 

High  Intensity  Residential 

0.5 

23 

Commercial/Industrial/Transportation 

0.5 

31 

Bare  Rock/Sand/Clay 

1 

32 

Quarries/Strip  Mines/Gravel  Pits 

1 

33 

Transitional 

1 

41 

Deciduous  Forest 

1 

42 

Evergreen  Forest 

1 

43 

Mixed  Forest 

1 

51 

Shrubland 

1 

61 

Orchards/Vineyards/Other 

1 

71 

Grasslands/Herbaceous 

1 

81 

Pastures/Hay 

1 

82 

Row  Crops 

1 

83 

Small  Grains 

1 

84 

Fallow 

1 

85 

Urban/Recreational  Grasses 

0 

91 

Woody  Wetlands 

0 

92 

Emergent  Herbaceous  Wetlands 

0 
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3.3.19  intersection 

This  map  is  based  on: 

interstates  (page  17) 
otherroads  (page  17). 

The  steps  of  the  algorithm  consist  of: 

1.  Find  all  roads  categorized  as  county  or  state. 

2.  Thin  the  results  to  make  sure  that  no  road  is  wider  than  a  single 
pixel. 

3.  Find  all  road  pixels  in  the  results  that  border  more  than  two  other 
road  pixels  (i.e.,  intersections),  and  save  the  results. 

3.3.20  hwyBuff 

This  map  is  based  on: 

interstates  (page  17) 
otherroads  (page  17). 

The  purpose  of  this  analysis  is  to  ensure  no  traffic  movement  off  limited- 
access  highways  except  at  intersections  with  other  roads.  (Note  that  this 
analysis  is  associated  with  an  unresolved  challenge  that  involves  restrict¬ 
ing  access  to  these  highways  -  except  at  ramps.  That  is,  traffic  should  be 
allowed  to  flow  under/over  limited  access  highways  without  connecting  to 
them. 

The  steps  of  the  algorithm  consist  of: 

1.  All  road  cells  are  given  category  o  (no  buffer). 

2.  All  cells  that  border  road  category  1  (limited  access  highway)  cells 
are  given  category  1  (buffer). 

3.3.21  developable 

This  map  is  based  on: 

landcover  (page  16) 
noGrowth_att  (page  16). 

This  is  simply  a  convenient  map  that  is  used  by  a  number  of  other  analy¬ 
ses,  which  indicates  whether  or  not  cells  can  be  developed  for  residential. 
Cells  marked  as  no-growth  or  containing  landcover  categories  11,  23,  91, 
92,  or  85  are  given  category  o  (no-development).  Cells  containing  land- 
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cover  categories  21,  22,  and  23  are  given  category  2  (developed),  and  the 
remaining  cells  are  developable  (category  1). 

3.4  Development  of  Probabilities  Map 

The  above  attractor  maps  mostly  contain  factual  information  that  must 
now  be  converted  to  development  probability  maps.  The  algorithm  for 
each  is  identical  and  provides  a  probability  of  finding  residential  areas 
within  each  category  in  each  attractor  map.  The  result  is  a  probability 
graph  relating  attractor  map  category  values  to  probability  of  finding  resi¬ 
dential  areas.  This  graph  is  then  applied  to  the  attractor  map  to  create  a 
probability  map. 

The  implemented  algorithm  works  as  follows: 

1.  The  attractor  map  is  normalized. 

a.  The  highest  category  in  the  map  is  determined. 

b.  The  map  is  divided  by  this  highest  value  and  multiplied  by  10.  The 
result  is  an  integer  map  with  values  of  o  to  10  that  represent  ranges 
of  the  original  categories. 

2.  This  result  is  cross-tabulated  with  the  “developed”  map,  which  results 
in  a  coincident  tabulation  of  the  categories  in  each  map. 

3.  For  each  of  the  10  categories,  the  number  of  developed  areas  are  di¬ 
vided  by  the  sum  of  the  number  of  developed  and  developable  areas  to 
give  the  probability  of  finding  developed  areas  in  each  category,  which 
produces  a  graph.  (Figure  9  shows  an  example.) 

4.  The  resulting  graph  is  applied  to  the  attractor  map  to  create  the  prob¬ 
ability  map.  The  straight  line  overlaid  on  the  histogram  in  the  figure 
represents  the  graph  applied. 
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Figure  9.  Sample  graph  for  converting  an  attractor  to  a  probability  of  urban  area. 

The  end  result  is  a  factual  map  showing  the  probability  of  finding  urban 
areas  within  the  categories  of  the  attractor  map.  The  collection  of  probabil¬ 
ity  maps  (Table  1)  can  then  be  combined  to  generate  an  overall  probability 
map. 


Table  1.  Collection  of  probability  maps. 


Probability  Map  Name 

Based  on 

cities_prob 

cities_att  (page  17),  developable  (page  25) 

roacLprob 

road_att  (page  20),  developable  (page25) 

ramp_prob 

ramp_att  (page  18),  developable  (page  25) 

intersect_prob 

intersect_att  (pagel9),  developable  (page  25) 

tree_prob 

forest_att  (page20),  developable  (page  25) 

water_prob 

water_att  (page20),  developable  (page  25) 

slope_prob 

slope_att  (page  18),  developable  (page  25) 

staterd_prob 

staterd_att  (page  19),  developable  (page  25) 

3.5  Development  of  Final  Maps 

Three  useful  maps  can  now  be  generated  in  sequence.  The  first  is  a  resi¬ 
dential  attractor  map,  which  is  a  combination  of  the  probability  maps.  The 
second  is  a  probability  map  based  on  this  attractor  map  and  is  developed 
with  the  same  algorithm  described  above.  Finally,  using  this  map,  a  future 
pattern  of  land  development  can  be  generated. 
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3.5.1  residential  att 


This  map  is  based  on: 
developable  (page  25) 
slope_prob  (page  27) 
staterd_att  (page  27) 
water_prob  (page  27) 
intersect_prob  (page  27) 
landcover  (page  16) 


road_prob  (page27) 
ramp_prob  (page  27) 
tree_prob  (page  27) 
cities_prob  (page  27) 
noGrowth_att  (page  16) 
boundary  (page  16) 


This  provides  the  first  combination  of  all  of  the  selected  attractors  to  resi¬ 
dential  development.  The  values  are  simply  summed  and  divided  by  the 
total  number  of  values  to  give  a  result  in  the  range  of  0-1.  While  it  is  possi¬ 
ble  to  weight  these,  the  importance  of  each  is  fully  captured  in  the  process 
generating  each.  The  total  number  of  residential  areas  in  a  study  area  is 
represented  by  the  area  under  the  curve  in  Figure  9.  Regardless  of  the 
shape  of  the  curve  (see  samples  in  Figure  7,  p  13),  the  area  under  the  curve 
will  be  the  same.  Attractors  with  little  attraction  will  be  associated  with 
relatively  flat  graphs  and  will  therefore  have  little  consequence  on  the  re¬ 
sults  of  combining  the  various  probabilities  associated  with  all  attractors. 
Therefore,  different  weights  are  not  used,  nor  are  any  justified  with  this 
simple  approach. 


3.5.2  residential_prob 

This  map  is  based  on: 

residential_att  (page27) 
developable  (page  25). 


The  above  residential  attractor  map  (residential_att)  can  now  be  proc¬ 
essed  with  the  identical  algorithm  associated  with  Figure  9,  above.  This  is 
a  final  product  that  provides  the  best  available  picture  of  the  attractiveness 
to  urban  growth  across  the  study  area  based  on  the  statistical  evaluation  of 
the  attractiveness  of  chosen  factors. 

3.6  Assumptions  and  Caveats 

The  algorithm  described  has  a  number  of  important  assumptions  and  ca¬ 
veats.  It  was  developed  to  address  a  specific  question  within  time  and 
budget  constraints.  Each  important  assumption  is  discussed  below. 
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Only  Certain  User  Questions  Are  Answered 

LEAMram™  was  developed  to  identify  the  attractiveness  of  areas  within 
counties  to  future  urban  residential  development.  It  is  designed  to  help 
regional  planners  understand  regional  growth  potentials.  Other  ap¬ 
proaches  and  models  must  be  employed  to  predict  details  of  urban  growth 
within  cities  or  to  project  future  urban  growth. 

The  Current  Urban  Pattern  Is  Adequate  for  Predicting  Future  Urban 
Residential  Attractiveness 

Towns  and  cities  develop  over  very  long  periods  of  time  in  response  to 
many  technical,  economic,  social,  and  geographic  drivers.  Resulting  pat¬ 
terns  tend  to  persist,  but  new  development  can  be  different  due  to  techno¬ 
logical  developments  -  especially  those  in  transportation  and  communica¬ 
tion.  LEAMram™  analyzes  the  current  land  use  patterns  -  assuming  that 
they  represent  the  current  and  future  settlement  preferences.  Models  such 
as  LEAM™  analyze  land  use  change  or  registered  housing  starts  as  a  better 
representation  of  current  preferences.  Challenges  remain  to  predict  the 
future  residential  development  preferences  in  response  to  future  transpor¬ 
tation  costs  and  communication  potentials  offered  through  significantly 
enhanced  internet  communication  opportunities. 

The  Input  Map  Files  Are  Accurate 

Like  all  computer  software  analyses,  the  results  are  based  on  a  presump¬ 
tion  of  the  accuracy  of  the  inputs.  This  process  relies  on  nationally  avail¬ 
able  data  sets.  In  particular,  the  National  Land  Cover  Data  set  provides 
land  cover  classifications  at  a  resolution  of  30  meters  across  the  United 
States.  This  centrally  developed  data  set  is  primarily  based  on  image  proc¬ 
essing  and  may,  for  some  applications,  have  an  inadequate  level  of  accu¬ 
racy.  Of  course,  local,  better  data  may  be  used. 

Road  Systems  Cannot  Be  Overloaded 

A  most  important  part  of  the  LEAMram™  analysis  is  a  set  of  predictions  of 
the  travel  time  from  every  30-meter  square  parcel  to  important  attractors 
(intersections,  cities,  highways,  water,  etc.).  Roads  are  associated  with 
fixed  travel  times,  thereby  making  the  assumption  that  road  capacity  will 
not  be  reached,  or  if  reached  will  be  increased  through  a  construction  re¬ 
sponse. 
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There  Are  No  Differences  among  Owners  and  Parcels 

No  characteristics  of  owners  are  considered  including  wealth,  age,  phi¬ 
losophies,  primary  business,  etc.  In  reality,  motivations  and  philosophies 
of  owners  can  vary  significantly  resulting  in  different  land  use  desires. 
Similarly,  differences  in  parcels,  such  as  size,  are  not  considered.  Devel¬ 
opment  may  be  more  likely  on  large  parcels  that  can  be  purchased  through 
negotiations  with  a  single  owner  rather  than  on  an  area  that  involves  many 
parcels  and  owners.  Residential  development  is  more  cost  effective  for  lar¬ 
ger  developments. 

The  Value  of  Land  for  Other  Purposes  Is  Fixed 

This  is  an  important  assumption  that  will  not  be  true.  The  value  of  land  is 
the  amount  that  the  highest  bidder  is  willing  to  pay.  Land  may  convert  to 
residential  when  the  value  of  the  land  for  residential  (and  business/ind¬ 
ustrial)  rises  above  its  value  in  forest  or  agriculture.  However,  it  may  con¬ 
vert  if  the  latter  value  drops.  For  example,  a  region  that  traditionally  grows 
a  crop  that  drops  in  value  may  see  more  development  than  a  region  with 
solid  or  rising  crop  prices. 

Development  Attractiveness  Is  Inherent  in  Existing  Urban  Patterns 

That  is,  where  people  want  to  live  in  the  future  is  reflected  in  the  settle¬ 
ment  patterns  of  the  past.  Perhaps  a  more  accurate,  but  more  expensive 
approach  is  to  not  make  this  assumption.  Instead  of  calibrating  the  model 
on  the  complete  development  patterns,  calibrate  on  only  the  recent  devel¬ 
opment  derived  from  construction  start  data  or  on  differences  between 
consecutive  land  use  maps  -  separated  by  a  decade  or  more.  Only  the 
study  area  is  taken  into  account.  There  are  attractors  outside  the  study 
area  that  can  affect  the  study  area.  However,  despite  these  assumptions 
and  caveats,  LEAMram™  results  can  be  very  useful  in  a  regional  context  to 
identify  the  relative  attractiveness  of  development  in  a  rapidly  developing 


area. 
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4  The  Approach  Implementation 

The  algorithms  described  in  Chapter  2  have  been  implemented  using  the 
GRASS  GIS,  ver.  5.3.  GRASS  was  originally  developed  though  a  collabora¬ 
tive  effort  by  various  U.S.  government  agencies,  universities,  and  private 
companies  (Westervelt,  et  al.,  1987).  Active  GRASS  development  currently 
involves  universities,  agencies,  and  corporations  across  the  world  (Netler 
and  Mitasova,  2004).  GRASS  remains  a  freely  distributable  software  pack¬ 
age  that  runs  under  UNIX  and  Linux  systems.  GRASS  is  particularly  suit¬ 
able  for  the  cost-effective  support  of  the  LEAMram™  project  because  of 
the  wealth  of  raster  GIS  analysis  routines,  its  freely  distributable  nature, 
and  the  ability  to  write  any  necessary  new  software  by  adapting  existing 
software  source  code. 

The  primary  implementation  of  the  algorithms  is  through  a  UNIX  “make¬ 
file.”  The  “make”  program  reads  file  dependency  information  from  a  file, 
the  makefile,  and  executes  instructions  necessary  to  create  missing  or  out 
of  data  files  by  running  specific  UNIX  instructions  on  the  dependency 
files.  For  example,  consider  the  following  makefile: 

target:  file-a  file-b 
cat  file-a  file-b  >  target 
file-a: 

echo  “hello  ”  >  file-a 
file-b: 

echo  “world  ”  >  file-b 

Invoking  the  make  program  on  this  makefile  causes  the  program  to  iden¬ 
tify  the  existence  and  date  of  creation  of  target,  file-a,  and  file-b.  If  file-a 
and  file-b  exist,  and  the  creation  date  on  “target”  is  later  than  file-a  and 
file-b,  nothing  happens.  If  the  creation  date  is  later,  or  the  “target”  does 
not  exist,  the  following  command  is  executed:  “cat  file-a  file-b  >  target,” 
creating  “target.”  If  file-a  or  file-b  do  not  exist,  the  instructions  following 
their  entries,  respectively,  are  executed.  These  have  no  dependencies;  only 
their  existences  are  checked.  In  summary,  invoking  the  make  program  on 
this  file  will  result  in  three  files  being  created,  in  this  case  with  text. 

The  algorithms  described  in  the  previous  chapter  were  captured  in  a 
makefile  reproduced  in  Appendix  A.  The  general  format  of  this  file  is  based 
on  the  simple  example  above.  Comments  are  preceded  on  lines  in  the 
makefile  with  a  pound  (#)  sign.  You  are  encouraged  to  read  through  this 
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file,  cross-referencing  with  the  previous  algorithm  descriptions.  Note  that 
the  GRASS  GIS  programs  are  command-line  UNIX  programs  that  mix 
powerfully  with  standard  UNIX  commands  such  as  cat,  grep,  sed,  awk, 
and  others. 

4.1  Prerequisites 

The  following  prerequisites  must  be  addressed  to  run  the  LEAMram™ 
analysis.  First,  LEAMram™  is  GRASS  based  and  GRASS  is  primarily  avail¬ 
able  for  UNIX  and  Linux  systems.  Versions  of  GRASS  are  available 
through  the  Internet  free  of  charge.  GRASS  5.3  or  better  is  required  for 
this  application.  Acquire  and  install  GRASS  on  a  Linux  computer  of  your 
choice  and  become  familiar  with  Linux  and  GRASS  commands.  Also,  the 
make  program  must  be  available  on  your  Linux  system.  It  is  often  part  of 
Linux  software  development  packages,  which  may  not  have  been  initially 
installed. 

4.2  Create  Maps 

Acquire  the  basic  nationally  available  maps  from  the  Internet.  At  this 
point,  any  raster-based  GIS  of  your  choice  can  be  employed  to  import  the 
maps  into  a  common  coordinate  system  and  projection.  First,  create  the 
LEAMram™  base  input  maps,  adhering  to  the  map  names  and  contents. 
Then,  bring  the  maps  into  a  GRASS  database.  Geographic  locations  in 
GRASS  are  stored  as  sets  of  maps  called  mapsets.  Locations  are  directories 
(folders)  in  the  Linux  operating  system  and  mapsets  are  folders  within  the 
location.  If  needed,  create  a  new  GRASS  location,  which  will  create  the 
“PERMANENT”  mapset.  Ensure  that  the  projection  and  coordinate  system 
matches  that  chosen  when  the  maps  were  input  and  processed.  Once  a  lo¬ 
cation  and  mapset  are  created,  run  GRASS,  selecting  them.  Import  your 
base  maps  into  GRASS  using  any  of  a  variety  of  GRASS  map  import  com¬ 
mands,  depending  on  the  format  of  your  maps.  You  may  need  to  output 
your  maps  from  the  GIS  used  to  initially  process  the  maps  into  a  format 
“understood”  by  one  of  the  GRASS  import  commands. 

4.3  Run  LEAMram™ 

Create  a  makefile  with  the  contents  of  Appendix  A.  (We  recommend  that 
you  contact*  the  author  of  this  report  for  a  latest  version  of  this  makefile  as 
the  file  was  under  development  at  the  writing  of  this  document.)  Call  the 


*  Contact  information  is  available  by  searching  on  the  title  of  this  report  through  URL: 
http  ://www.  ce  ce  r.  a  rm  v.  m  i  l/te  c  h  re  po  rts 
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file  “makefile”  and  place  it  into  the  “cell”  directory  of  the  current  mapset. 
You  can  identify  the  directory  location  of  the  cell  file  by  looking  at  the 
GRASS  environment  variables  returned  by  the  g.gisenv  program.  Change 
directories  to  the  database,  then  the  location,  and  finally  the  mapset.  Then 
change  to  the  “cell”  directory  at  that  level.  Place  the  makefile  in  this  direc¬ 
tory. 

Edit  the  makefile  and  locate  the  names  of  the  base  maps  under  the  com¬ 
ment  “Provide  names  of  the  input  maps.”  Using  GRASS  display  com¬ 
mands,  review  these  maps  and  check  to  make  sure  they  have  contents  that 
match  their  descriptions  at  the  top  of  the  makefile.  Once  satisfied,  run  the 
command  “make”  while  in  the  “cell”  directory.  You  will  get  a  response 
similar  to  the  following: 

make  LEmaps  -  to  create  LUC  input  maps 

make  indexmaps  -  to  create  LUCs  index  maps 

make  residential_att  -  to  create  LUCs  residential  attractor  map 

make  residential_prob  -  to  create  LUCs  residential  probability  map 

make  futureLanduse  -  to  create  a  future  landuse  map 

make  clean  -  to  remove  temporary  maps 

make  veryclean  -  to  remove  all  maps  generated  by  this  makefile 


These  are  some  of  the  major  targets.  If  all  goes  well,  you  will  be  able  to 
successfully  run  the  command:  “make  residential_prob.”  The  entire  proc¬ 
ess  can  take  many  hours  depending  on  the  speed  of  your  computer,  its 
amount  of  main  memory,  and  the  size  of  your  study  area. 


ERDC/CERL  TR-06-28 


34 


5  A  Sample  Application 

Fort  Bragg,  NC  and  its  surrounding  communities  and  counties  are  eco¬ 
nomically  and  environmentally  linked.  Anticipated  changes  in  mission, 
weapon  systems,  and  training  will  require  substantially  larger  training  ar¬ 
eas  at  a  time  when  anticipated  population  growth  will  result  in  increased 
demands  for  urban  land.  These  demands  require  closer  analysis  of  the 
landscape,  predictions  of  future  needs,  identification  of  potential  conflicts, 
and  analysis  of  the  direct  and  indirect  consequences  of  proposed  local  in¬ 
vestments  and  policies.  The  effort  described  here  is  an  approach  for  identi¬ 
fying  areas  associated  with  a  significant  probability  of  urban  growth  in  the 
next  decades. 

Fort  Bragg,  like  many  military  installations  with  training  and  testing  mis¬ 
sions,  was  placed  in  a  relatively  remote  area  where  few  or  no  human  set¬ 
tlements  nearby  would  be  adversely  affected  by  the  military  activities,  and 
where  the  large  contiguous  tracts  of  land  necessary  to  support  assorted 
large-scale  activities  were  readily  and  inexpensively  accessible.  Purchased 
property  boundaries  were  demarcated  with  fences  and  signs.  Military 
training  and  testing  could  then  take  place  “within  the  fence  line.”  How¬ 
ever,  some  impacts  of  those  activities  extended  beyond  the  fence  line  in¬ 
cluding  noise,  dust,  smoke,  and  radio  transmissions.  Also,  in  the  interven¬ 
ing  years  human  settlements  and  agriculture  have  impacted  important 
habitats  and  populations  of  threatened  or  endangered  species  -  leaving 
the  installation  with  an  increasing  percentage  of  the  overall  habitat  and 
populations.  This  has  resulted  in  growing  concerns  about  community  and 
public  interest  in  the  potential  reduction  of  on-base  training  and  testing. 

The  procedures  described  in  the  previous  chapter  were  applied  to  the  Fort 
Bragg  area,  which  included  all  local  counties  participating  in  the  Sustain¬ 
able  Sandhills  effort.  The  following  images  show  a  study  area  consisting  of 
nearly  to  million  30-meter  square  grid  cells  (parcels)  across  the  multi¬ 
county  area.  Fort  Bragg  is  centered  in  the  middle  of  the  images.  The  pro¬ 
jection  is  Albers  Equal  Area  resulting  in  the  north  direction  being  slightly 
to  the  left  of  vertical.  Each  image  is  associated  with  a  color  table  spanning 
the  scale  from  o  to  1.  Note  that  some  images  are  richer  in  color  indicating  a 
tighter  correlation  between  urban  areas  and  the  particular  driver. 
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The  production  of  the  “attraction  to  cities”  map  (Figure  to)  is  significantly 
different  from  the  creation  of  the  other  maps  (Figures  n  through  16).  The 
goal  is  to  produce  a  map  that  captures  the  attractiveness  to  jobs,  shopping, 
and  friends  as  a  function  of  the  innate  attractiveness  of  those  qualities  me¬ 
diated  by  distance  from  them.  Production  of  this  map  in  an  automated  way 
that  can  be  readily  applied  to  any  location  is  challenging.  The  first  chal¬ 
lenge  is  to  identify  the  areas  or  points  that  represent  centers  of  attraction, 
which  is  accomplished  during  the  manual  preparation  process.  Ideally  a 
separate  analysis  would  be  done  for  each  small  town  and  for  each  city’s 
neighborhoods  and  employment  centers.  This  has  been  determined  to  be 
impractical.  The  current  approach  is  to  divide  cities  into  four  categories: 
small,  medium,  large  and  extra-large.  Each  is  associated  with  an  average 
population  for  the  group,  which  provides  the  strength  of  attraction.  Cen¬ 
ters  of  this  attraction  are  then  determined  for  each  separate  site  and  are 
typically  the  densest  area  near  the  center  of  the  city  or  town.  The  auto¬ 
matic  processing  then  consists  of  computing  minimal  travel  times  from 
each  cell  in  the  analysis  area  to  the  nearest  center  of  attraction  in  each 
category.  These  are  converted  to  attractiveness  values,  which  are  the  popu¬ 
lation  size  divided  by  the  square  of  the  distance  from  the  attraction  center. 


Figure  10.  Attraction  to  cities. 
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Figure  11.  Attraction  to  roads. 


Figure  12.  Attraction  to  ramps  on  limited  access  highways. 
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Figure  13.  Attraction  to  major  intersections. 


Figure  14.  Attraction  to  different  slopes. 
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Figure  16.  Attraction  to  lots  bordering  water. 
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The  combination  of  the  attractors,  displayed  in  Figure  17,  show  significant 
attraction  to  new  growth  around  Fayetteville  and  the  Dunn/Erwin  area  lo¬ 
cated  in  the  northeast  corner.  Zooming  into  the  southern  boundary  of  the 
installation  and  the  western  edge  of  Fayetteville,  reveals  significant  prefer¬ 
ence  areas  close  to  the  city.  However,  some  attraction  is  visible  associated 
with  existing  rural  neighborhoods  closer  to  the  city. 

The  final  map  created  in  this  effort  (Figure  18)  provides  a  relative  index 
value  for  every  location  representing  its  attractiveness  to  residential 
growth  after  considering  a  number  of  location  factors  listed  again  here: 

•  location  within  the  study  area 

•  whether  residential  growth  is  possible 

•  current  land  use 

•  distance  to  water 

•  distance  to  forest 

•  slope 

•  distance  to  cities  (city  centers  /  employment  centers) 

•  distance  to  county  road 

•  distance  to  state/federal  highway 

•  distance  to  limited  access  highway 

•  distance  to  nearest  intersection. 


Figure  17.  Summary  attraction  to  residential  development. 
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Figure  18.  Summary  attraction  along  southern  border  of  Fort  Bragg. 

Several  assumptions  and  caveats  need  to  be  recognized: 

•  Other  potential  attraction  considerations  could  be  added  to  this  list  and 
some  on  the  list  could  be  dropped  for  any  particular  location.  Note,  for 
example,  that  drinking  water  availability  is  not  considered  and  can  be 
critical  from  a  legal  standpoint  (e.g.,  western  water  rights)  or  from  a 
geologic  perspective.  Some  areas  offer  more  well  water  opportunities 
than  others. 

•  Zoning  is  not  considered  in  this  analysis  and  can  be  important  in  the 
attraction  of  urban  growth. 

•  This  analysis  assumes  no  new  investments  in  roads  or  in  the  develop¬ 
ment  of  new  neighborhoods,  but  can  be  rerun  with  such  developments 
provided  as  inputs.  It  also  does  not  recognize  the  affect  of  development 
affecting  travel  times  on  roads. 

•  Travel  times  are  assumed  to  be  optimal  for  each  type  of  road.  Hence, 
the  resulting  maps  provide  a  snapshot  in  time  of  the  attraction  to  new 
growth,  but  do  not  consider  the  impact  of  new  growth  on  the  overall  at¬ 
traction. 

•  The  analysis  identifies  attractiveness  to  new  development  on  a  cell-by- 
cell  basis  -  parcels  approximately  the  size  of  city  lots.  While  some 
growth  happens  this  way,  much  development  occurs  as  part  of  new 
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neighborhood  development  sites  that  can  be  roughly  800  meters 
square. 

•  Parcel  size  and  ownership  is  not  considered.  Developers  looking  to 
build  a  new  neighborhood  are  more  likely  to  purchase  a  single  large 
parcel  rather  than  piece  together  many  smaller  contiguous  parcels  - 
making  urban  development  less  likely. 

•  Negative  attractors  are  not  considered  in  this  analysis.  Urban  devel¬ 
opment  tends  to  avoid  being  co-located  with  industrial  sites. 

•  The  attractiveness  to  urban  growth  can  be  outbid  by  attractiveness  to 
other  land  uses  such  as  parks  and  industrial  areas.  This  competition  is 
not  identified  here. 

•  The  attractiveness  to  cities  map  must  be  developed  in  close  consulta¬ 
tion  with  local  planners  to  best  capture  understandings  of  the  location 
and  attractiveness  of  local  population,  employment,  and  shopping  cen¬ 
ters. 
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6  Summary  and  Recommendations 

This  work  has  demonstrated  an  approach  (process  and  procedure)  that 
can  be  used  in  close  coordination  with  regional  planners  to  test  the  conse¬ 
quences  of  alternative  regional  plans  using  an  inexpensive  GIS-based  ap¬ 
proach  for  predicting  future  urban  growth  patterns  that: 

1.  uses  nationally  available  data 

2.  requires  minimal  data  preparation 

3.  is  readily  calibrated 

4.  does  not  require  time-series  data 

5.  identifies  the  relative  attractiveness  of  areas  around  military  installa¬ 
tions  much  better  than  the  current  image-processing  based  trend 
analyses  and  much  less  expensively  than  currently  available  urban 
growth  simulation  models. 

(Note  that  results  generated  by  the  approach  are  not  ready  to  be  used  in 
informing  local  planning  processes.) 

This  effort  is  part  of  a  larger  effort  that  includes  the  development  of  a 
simulation  model  to  project  future  landscape  settlement  patterns.  The  lar¬ 
ger  effort  is  called  the  Land  use  Evolution  and  impact  Assessment  Model 
(LEAM™)  and  supports  the  development  of  the  LEAM  Land  Use  Change 
(LEAM-LUC™)  model.  This  model  and  the  larger  effort  are  addressing 
some  of  the  caveats  associated  with  the  limited  analysis  presented  here.  In 
particular,  LEAM-LUC™  considers  the  effect  of  development  on  future  de¬ 
velopment  and  considers  competition  for  turning  undeveloped  areas  into 
residential,  commercial,  and  open  space. 

The  next  available  steps  are: 

1.  Review  the  analysis  process  with  local  planners  and  modify  inputs  to 
better  capture  the  local  planning  scenarios.  This  will  allow  creation  of  a 
development  attractiveness  map  that  best  matches  the  local  urbaniza¬ 
tion  dynamics. 

2.  Test  alternative  planning  scenarios  that  might  include  upgrading  roads 
and  highways  or  building  new,  closing  roads,  and  zoning  ideas. 

3.  Run  the  LEAM-LUC™  model  to  generate  future  scenarios  that  could 
develop  in  direct  and  indirect  response  to  economic,  population,  and 
employment  conditions. 
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Appendix  A:  Unix  Input  File  To  Create  LEAM 
Landuse  Evolution  Input  Files  Using 
GRASS  Commands 

############################################################################# 

#  last  modified:  Tue  Apr  20  11:25:57  CDT  2004 

# 

#  This  file  is  used  by  the  unix  make  program  to  create  the 

#  LEAM  landuse  evolution  input  files  using  grass  commands. 

# 

#  To  use, 

#  Start  GRASS  and  use  a  mapset  that  will  contain  the  results 

#  Create  a  directory  named  cell  in  this  mapset  (if  not  already  there) 

#  Copy  this  Makefile  to  that  directory 

#  Set  the  name  of  the  needed  input  maps  (see  below) 

#  Set  the  size  of  the  cities  (see  below) 

#  Set  the  attractor  weights  (see  below) 

#  From  the  cell  directory,  run  make  on  this  file 

############################################################################# 

#  Input  Maps  and  Map  Categories 

#  Make  sure  the  input  maps  use  the  following  categories  . . . 

#  SMALL_CITY  -  Location  of  "small  cities" 

#  0  No  small  city 

#  1  Small  city 

#  MED_CITY  -  Location  of  "medium  cities" 

#  0  No  small  city 

#  1  Small  city 

#  LARGE_CITY  -  Location  of  "large  cities" 

#  0  No  small  city 

#  1  Small  city 

#  NO_GRCWTH  -  Explicit  identification  of  areas  of  no  growth 

#  0  No  restrictions 

#  1  No  growth  areas 

#  DEM 

#  Categories  are  meters  above  sea  level 

#  M(JNICIPAL_BOUNDARY 

#  City  boundaries 

#  STUDY_AREA 

#  0  Outside  study  area 

#  1  Inside  study  area 

#  INTERSTATES 

#  0  No  road 

#  1  Limited  access  highway 

#  CTHERROADS 

#  0  No  road 

#  2  Federal  highway 

#  3  State  highway 

#  4  County/neighborhood  road 

#  5  (New?)  county/ neighborhood  road 

#  6  Ramps 

#  7  Private  road 

#  8  Neighborhood  road 

#  LANDCOVER  -  NLCD  Landcover  Map 

#  11  Open  water 

#  12  Perennial  Ice/Snow 

#  21  Low  Intensity  Residential 

#  22  High  Intensity  Residential 

#  23  Commercial/ Industrial/Transportation 
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#  31  Bare  Bock/Sand/Clay 

#  32  Quarries/Strip  Mines/Gravel  Pits 

#  33  Transitional 

#  41  Deciduous  Forest 

#  42  Evergreen  Forest 

#  43  Mixed  Forest 

#  51  Shrubland 

#  61  Orchards/Vineyards/Other 

#  71  Grasslands/Herbaceous 

#  81  Pastures/Hay 

#  82  Bow  Crops 

#  83  Small  Grains 

#  84  Fallow 

#  85  Urban/Becreational  Grasses 

#  91  Woody  Wetlands 

#  92  Emergent  Herbaceous  Wetlands 

# 

############################################################################# 

#  Provide  names  of  the  input  maps  (name@mapset) 

SMALL_CITY=small_city@leamBase 

MEDjGITY=medium_city@leamBase 

LABGE_CITY=large_city@leamBase 

XL_CITY=xl_city@leamBase 

NO_GBCWTH=nogrowth@learriBase 

DEM=den@learriBase 

MUNICIPAL_BOUNDARY=mnbnd@leamBase 

STUDY_ABEA=county@leamBase 

LANDCOVEB=landcover@leamBase 

INTEBSTATES=interstates@leamBase 

OTHEBROADS=otherroads@leamBase 

# 

############################################################################# 

#  Set  the  population  average  for  large,  medium,  and  small  cities.  These 

#  will  be  used  to  create  the  cities_attractor  map,  which  combines  these 

#  populations  driving  time  maps. 

XLPOP=150000 

LABGEPOP=95000 

MEDIUMPOP=50000 

SMALLPOP=7000 

POP2 000=75 654 4 
DENSITY=1 . 68 

# 

############################################################################# 

#  Set  the  multipliers  for  residential  growth  drivers.  These  are  multiplied 

#  times  their  respective  0-1  index  maps.  The  results  are  summed  and  divided 

#  by  the  sum  of  the  multipliers  to  give  an  overal  0-1  growth  attractor 

#  index  map. 

#  NOTE:  DOES  NOT  MATCH  THE  GEBTNEB,  ET  AL.  AND  LUC  EQUATION 
intersect_prob_MULT=l 

road_prob_MULT=l 

slope_prob_MULT=l 

ramp_prob_MULT=l 

staterd_prob_MULT=l 

forest_prob_MULT=l 

water_prob_MULT=l 

cities_prob_MULT=l 

neighbor_prob_MULT=l 

############################################################################# 

############################################################################# 

############################################################################# 

#  IDEALLY  NO  EDITING  IS  BEQUIBED  BELOW  HEBE 
############################################################################# 
############################################################################# 
############################################################################# 
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LEmaps=landcover  water_att  forest_att  road_att  ramp_att  cityBuff  intersect_att  staterd_att 
slope_att  noGrowth  boundary  cities_att2  neighbor_att  cities_att2 

IndexMaps=staterd_prob  slope_prob  water_prob  forest_prob  intersect_prob  \ 
rarrp_prob  road_prob  cities_prob 

CtherMaps=hwyBuff  road  travelTime30  landcover  municipalboundary  \ 

intersection  smallCity  mediurrCity  largeCity  xLargeCity  dan  travelSpeed30  travelTime90  \ 
xlargeCityTime  largeCityTime  mediurrCityiime  smallCityTime  \ 
landTravelSpeed30  \ 

travelTime  developable  residential_att  futureLanduse 

LEfiles=$ (shell  echo  ${tEmaps}  |  sed  -e  's/  /, /g') 

Otherfiles=$ (shell  echo  ${OtherMaps}  |  sed  -e  's/  /,/g') 

Indexfiles=$ (shell  echo  ${IndexMaps}  |  sed  -e  's/  /,/g') 

default : 

@echo  make  LErnaps  -  to  create  LUC  input  maps 
@echo  make  indexmaps  -  to  create  LUCs  index  maps 

@echo  make  residential_att  -  to  create  LUCs  residential  attractor  map 
@echo  make  residential_prob  -  to  create  LUCs  residential  probability  map 
@echo  make  futureLanduse  -  to  create  a  future  landuse  map 
@echo  make  clean  -  to  remove  temporary  maps 

@echo  make  veryclean  -  to  remove  all  maps  generated  by  this  makefile 
LQraps:  ${LBraps} 
indexmaps:  ${IndexMaps} 

###################################################################### 

###################################################################### 

#  Fetching  the  base  maps 

###################################################################### 

###################################################################### 

#  Fetches  a  map  of  the  locations  of  the  small  cities 
smallCity: 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  ' smallCityO=if (${SMALL_CITY}  <=  0,null(),  ${SMALL_CITY} ) '  >  /dev/null  2>&1 
r.clurrp  smallCityO  out=smallCityl 

#  -  Find  centroids 

r.mapcalc  ' smallCity2=if (isnull (smallCityl) ,null () , 1) '  >  /dev/null  2>&1 

r.  volume  data=smallCity2  clurrp=smallCityl  site=smallCityCentroids 

#  -  convert  to  cells 

s. to.rast  smallCityCentroids  out=smallCity3  size=3 
r.null  smallCity3  setnull=0 

r.mapcalc  ' smallCity=if (isnull (smallCity3) , null (), 1) '  >  /dev/null  2>&1 
g . remove  smallCityO , smallCityl , smallCity2 , smallCity3 

#  Fetches  a  map  of  the  locations  of  the  medium  cities 
mediurrCity: 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  'mediumCityO=if (${MED_CITY}  <  0,0,  ${MED_CITY} ) '  >  /dev/null  2>&1 
r.clurrp  mediumCityO  out=mediurrCityl 

#  -  Find  centroids 

r.mapcalc  'mediurrCity2=if  (isnull  (mediurrCityl) , null  (),  1)  '  >  /dev/null  2>&1 

r.  volume  data=mediumCity2  clurrp=mediumCityl  site=mnediurrCityCentroids 

#  -  convert  to  cells 

s. to.rast  mediumCityCentroids  out=mediumCity3  size=3 
r.null  mediumCity3  setnull=0 


ERDC/CERL  TR-06-28 


47 


r.mapcalc  'mediurrCity=if  (isnull  (mediumCity3) , null  (),  1)  '  >  /dev/null  2>&1 
g .  remove  mediurrCityO ,  mediurrCity  1 ,  mediurrCity2 ,  mediurrCity3 

#  Fetches  a  map  of  the  locations  of  the  large  cities 
largeCity: 

@echo  1 ######################### 1 ;  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  ' largeCityO=if ( $ { LfiRGE_CITY }  <  0,0,  ${LARGE_CITY} ) '  >  /dev/null  2>&1 
r.clurrp  largeCityO  out=largeCityl 

#  -  Find  centroids 

r.mapcalc  ' largeCity2=if (isnull (largeCityl) , null (), 1) '  >  /dev/null  2>&1 

r.  volume  data=largeCity2  clump=largeCityl  site=largeCityCentroids 

#  -  convert  to  cells 

s. to.rast  largeCityCentroids  out=largeCity3  size=3 
r.null  largeCity3  setnull=0 

r.mapcalc  'largeCity=if (isnull (largeCity3) , null (), 1) '  >  /dev/null  2>&1 
g . remove  largeCityO , largeCityl , largeCity2 , largeCity3 

#  Fetches  a  map  of  the  locations  of  the  xl  cities 
xLargeCity: 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  'xLargeCityO=if (${XL_CITY}  <  0,0,  ${XL_CITY}) '  >  /dev/null  2>&1 
r.clurrp  xLargeCityO  out=xIargeCityl 

#  -  Find  centroids 

r.mapcalc  'xLargeCity2=if (isnull (xlargeCityl) , null (), 1) '  >  /dev/null  2>&1 

r.  volume  data=xLargeCity2  clurrp=xLargeCityl  site=xIargeCityCentroids 

#  -  convert  to  cells 

s. to.rast  xLargeCityCentroids  out=xLargeCity3  size=3 
r.null  xLargeCity3  setnull=0 

r.mapcalc  'xLargeCity=if (isnull (xLargeCity3) ,null (), 1) '  >  /dev/null  2>&1 
g . remove  xLargeCityO , xlargeCityl , xIargeCity2 , xLargeCity3 

#  Fetches  a  map  of  no  growth  areas. 

#  Should  be  replaced  with  fetching  maps  of  various  no  growth  areas  that  can  be 

#  edited  and  then  merged. 
noGrowth: 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  'noGrowth=if (round (${NO_GROWTH}) ) '  >  /dev/null  2>&1 
r.null  noGrowth  null=0 

#  Fetches  the  digital  elevation  model 
dem: 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  'derrt=${DEM} '  >  /dev/null  2>&1 

#  Create  a  municipal  boundary  map  based  on  the  municipal  boundary  in  PERM 
rnunicipalboundary: 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r. buffer  in=${MUNICIPAL_BOUNDARY}  out=minicipalboundary  dist=1.5  unit=miles 

#  Fetch  the  study  area  boundary  frcm  PERM 
boundary: 

@echo  1 ######################### 1 ;  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  'boundary=if  (isnull  (${STUDY_AREA} )  ||  ${STUDY_AREA}  =0,null  () ,  1)  '  >  /dev/null 

2>&1 

#  Fetches  the  landcover  file  frcm  PERM  and  masks  it  with  the  study  area  boundary 
landcover:  boundary 
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@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

r .mapcalc  landcover= ' if (isnull (boundary) , null ( ) ,  \ 

if (isnull (${LSNDCOVER}) , 0, ${LANDCDVER} ) ) '  >  /dev/null  2>&1 
r.null  landcover  setnull=0 
(echo  11  110  127  177  ;  \ 
echo  21  253  229  228  ;  \ 
echo  22  247  178  159  ;  \ 
echo  23  229  86  78  ;  \ 
echo  31  210  205  192  ;  \ 
echo  32  175  175  177  ;  \ 
echo  33  83  62  118  ;  \ 
echo  41  133  199  126  ;  \ 
echo  42  56  129  78  ;  \ 
echo  43  212  231  176  ;  \ 

echo  51  220  202  143  ;  \ 

echo  61  187  174  118  ;  \ 

echo  71  253  233  170  ;  \ 

echo  81  251  246  93  ;  \ 

echo  82  202  145  70  ;  \ 

echo  83  121  108  74  ;  \ 

echo  84  244  238  202  ;  \ 

echo  85  240  156  54  ;  \ 

echo  91  200  230  248  ;  \ 

echo  92  100  179  213  ;  \ 

echo  end  ;  )  I  \ 
r. colors  landcover  color=rules 

#  Fetches  the  road  files 
otherroads:  boundary 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r. mapcalc  $@='if  (isnull  (boundary) , null  (), if  (${  OTHERROADS  }=0, null  (),  ${ OTHERROADS}) )  '  > 
/dev/null  2>&1 

r.null  $@  null=0 

interstates:  boundary 

@echo  1 ######################### 1 ;  echo  $@  ;  date 
g. region  res=30 

r  .mapcalc  $@='if  (isnull  (boundary) ,  null  (),  if  (${ INTERSTATES }==0,  null  (),  1) )  '  >  /dev/null 

2>&1 

r.null  $@  null=0 
cross:  otherroads  interstates 

r. mapcalc  cross='if  (interstates=l  &&  otherroads=6,  1,  0)' 

#  Force  state  highways  to  connect  to  interstates 
#road:  otherroads  interstates 

#  #  Identify  state  highway  cells  that  border  interstates  as  ramps 

#  r. mapcalc  'road=if (interstates>0, interstates, otherroads) ' 

#  r.null  road  null=0 

###################################################################### 

###################################################################### 

#  The  city  attractor  map  is  prepared  in  stages.  The  basic  approach  is 

#  as  follows: 

#  Fetch  the  small,  medium,  and  large  city  footprint  maps 

#  For  each  (small,  medium,  and  large) 

#  Create  road  travel  time  map  at  30  meter  resolution 

#  Create  overland  travel  time  map  starting  from  times  set  above 

#  Merge  the  two  maps 

#  Combine  the  small,  medium,  and  large  cumulative  maps 
###################################################################### 
###################################################################### 
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. ./LEAM/awkscript: 

@mkdir  -p  . . /LEAM 

@echo  'BEGIN  {  FS  =  "  cat=0;tot=0;developed=0;developable=0;gotone=0  }'  > 

. . /LEAM/awkscript 

@echo  ' { '  »  . . /LEAM/awkscript 

@echo  '  if  ($$1  !=  cat)'  »  ../LEAM/awkscript 

@echo  '  { '  »  . . /LEAM/awkscript 

@echo  '  if  (cat  !=  0  &&  cat  !=  "*") '  »  . . /LEAM/awkscript 
@echo  '  { '  »  . . /LEAM/awkscript 

@echo  '  if  (developable  !=  0) '  »  . . /LEAM/awkscript 

@echo  '  { '  »  . . /LEAM/awkscript 

@echo  '  if  (gotone  !=0)  printf (",  ")  ;'  »  ../LEAM/awkscript 

@echo  '  printf ("%f,%f",  cat*MULT/20,  developed/developable)  ;'  » 

. . /LEAM/awkscript 

@echo  '  gotone  +=!;'».. /LEAM/awkscript 

@echo  '  } '  »  . . /LEAM/awkscript 

@echo  '  } '  »  . . /LEAM/awkscript 

@echo  '  developed  =  0;  developable  =  0  ;  cat  =  $$1'  »  ../LEAM/awkscript 
@echo  '  } '  »  . . /LEAM/awkscript 

@echo  '  if  ($$2  =  2  | |  $$2  =  1)  developable  +=  $$3; '  »  ../LEAM/awkscript 

@echo  '  if  ($$2  =  2  )  developed  =  $$3; '  »  ../LEAM/awkscript 

@echo  '  tot  +=  $$3; '  »  ../LEAM/awkscript 

@echo  ' } '  »  . . /LEAM/awkscript 

@echo  'END  {printf ("\n") } '  »  ../LEAM/awkscript 

#  Combine  the  travel-time  maps  to  get  to  nearest  small,  medium,  and  large 

#  city  into  one  map  using  an  inverse-distance  weighting  algorithm  using  the 

#  population  of  each  type  as  the  base  weight. 

cities_att:  smallCityTime  mediumCityTime  largeCityTime  xLargeCityTime  developable 
@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  $@__tl=' round  (if  (developable,  \ 

float  (${XLPOP})  /if  (xLargeCityTime=0, 1, xLargeCityTime)  +  \ 
float  (${LARGEPOP}) /if  (largeCityiime=0,l,  largeCityTime)  +  \ 
float  ($  {MEDIUMPOP} )  /if  (mediumCityTime=0,  l,mediurrCityTime)  +  \ 
float  (${SMALLPOP} )  /if  (smallCityrime=0, 1,  smallCityTime) ,  0) )  '  \ 

>  /dev/ null  2>&1 

#  This  forces  a  4-byte  file  for  LEAMluc 

r.mapcalc  cities_att=' if  (row  ()=1  &&  col  0=1,33554432,  $@_tl)  '  \ 

>  /dev/ null  2>&1 
g. remove  $@_tl 

roadBuffer:  otherroads 

@echo  ' ######################### ' ;  echo  $@  ;  date 
r. buffer  -z  otherroads  out=$@  \ 

distance=30, 60, 90, 120, 150, 180, 210, 240, 270, 300, 330, 360, 390, 420, 450, 480, 510, 540, 570, 600,  630 

,660,690 

smallCityTimeRoad:  cross  intTravelTime30  othTravelTime30  smallCity 
@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

/data/FF/bin/r.cost.dks  in=othTravelTime30  m2=intTravelTime30  \ 
out=$@  xover=cross  start_rast=smallCity  percent=100  \ 

>  /dev/ null  2>&1 

#  Create  a  travel  cost  surface  for  travel  to  small  cities 
smallCityTime:  smallCityTimeRoad  overlandTravelTime90 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

r.null  smallCityTimeRoad  null=0 
r.mapcalc  '$@_tl  =  max(  \ 
smallCityTimeRoad [— 1,  —1] ,  \ 
smallCityTimeRoad [— 1, 0] , \ 
smallCityTimeRoad  [— 1, 1] , \ 
smallCityTimeRoad [  0 , -1 ] , \ 
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smallCityTimeRoad  , \ 
smallCit yTimeRoad [  0 , 1 ] , \ 
smallCityTimeRoad [  1, -1 ] , \ 
smallCityTimeRoad [  1 , 0 ] , \ 
smallCityTimeRoad [  1,1])'  >  /dev/null  2>&1 
r.null  $@_tl  setnull=0 
r.null  smallCityTimeRoad  setnull=0 
g. region  res=90 

r.cost  -r  start_rast=$@_tl  input=overlandTravelTime90  output=$@_t2  percent=100 
g. region  res=30 

r.mapcalc  smallCityTime='if (isnull (smallCityTimeRoad) , $@_t2, smallCityTimeRoad) ' 
g. remove  $@_tl,$@_t2 

mediumCityTimeRoad:  cross  intTravelTime30  othTravelTime30  mediumCity 
@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

/data/FF/bin/r.cost.dks  in=othTravelTime30  m2=intTravelTime30  \ 
out=$@  xover=cross  start_rast=mediumCity  percent=100  \ 

>  /dev/null  2>&1 

#  Create  a  travel  cost  surface  for  travel  to  medium  cities 
mediumCityTime:  mediumCityTimeRoad  overlandTravelTime90 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

r.null  mediumCityTimeRoad  null=0 
r.mapcalc  '$@_tl  =  max(  \ 
mediumCityTimeRoad [-1,-1] , \ 
mediumCityTimeRoad[-l, 0] , \ 
mediumCityTimeRoad [-1 , 1 ] ,  \ 
mediumCityTimeRoad [  0, -1] , \ 
mediumCityTimeRoad  , \ 
mediumCityTimeRoad [  0, 1 ] , \ 
mediumCityTimeRoad  [  1 , -1 ] , \ 
mediumCityTimeRoad  [  1 , 0 ] , \ 
mediumCityTimeRoad [  1,1])'  >  /dev/null  2>&1 
r.null  $@_tl  setnull=0 
r.null  mediumCityTimeRoad  setnull=0 
g. region  res=90 

r.cost  -r  start_rast=$@_tl  input=overlandTravelTime90  output=$@_t2  percent=100 
g. region  res=30 

r.mapcalc  mediumCityTime=' if (isnull (mediumCityTimeRoad) , $@_t2, mediumCityTimeRoad)  ' 
g. remove  $@_tl,$@_t2 

largeCityTimeRoad:  cross  intTravelTime30  othTravelTime30  largeCity 
@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

/data/FF/bin/r.cost.dks  in=othTravelTime30  m2=intTravelTime30  \ 
out=$@  xover=cross  start_rast=largeCity  percent=100  \ 

>  /dev/null  2>&1 

#  Create  a  travel  cost  surface  for  travel  to  large  cities 
largeCityTime :  largeCityTimeRoad  overlandTravelTime90 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

r.null  largeCityTimeRoad  null=0 
r.mapcalc  '$@_tl  =  max(  \ 
largeCityTimeRoad [-1, -1] , \ 
largeCityTimeRoad [-1 , 0] , \ 
largeCityTimeRoad [-1, 1] , \ 
largeCityTimeRoad [  0 , -1 ] , \ 
largeCityTimeRoad  , \ 
largeCityTimeRoad [  0 , 1 ] , \ 
largeCityTimeRoad [  1, -1] , \ 
largeCityTimeRoad [  1 , 0 ] , \ 
largeCityTimeRoad [  1,1])'  >  /dev/null  2>&1 
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r.null  $@_tl  setnull=0 

r.null  largeCityTimeRoad  setnull=0 

g. region  res=90 

r.cost  -r  start_rast=$@_tl  input=overlandTravelTime90  output=$@_t2  percent=100 
g. region  res=30 

r.mapcalc  largeCi tyTime=' if (isnull (largeCityTimeRoad) , $@_t2, largeCityTimeRoad) ' 
g. remove  $@_tl,$@_t2 

xlargeCityTimeRoad:  cross  intTravelTime30  othTravelTime30  xlargeCity 
@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

/data/FF/bin/r.cost.dks  in=othTravelTime30  m2=intTravelTime30  \ 
out=$@  xover=cross  start_rast=xLargeCity  percent=100  \ 

>  /dev/ null  2>&1 

#  Create  a  travel  cost  surface  for  travel  to  xlarge  cities 
xLargeCityTime:  xlargeCityTimeRoad  overlandTravelTime90 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.null  xlargeCityTimeRoad  null=0 
r.mapcalc  '$@_tl  =  max(  \ 
xLargeCityTimeRoad[-l,  -1] , \ 
xLargeCityTimeRoad[-l, 0] , \ 
xlargeCityTimeRoad [-1 , 1 ] ,  \ 
xLargeCityTimeRoad  [  0, —1] , \ 
xlargeCityTimeRoad  , \ 
xlargeCityTimeRoad [  0, 1 ] , \ 
xlargeCityTimeRoad [  1, -1] , \ 
xlargeCityTimeRoad  [  1 , 0 ] , \ 
xlargeCityTimeRoad [  1,1])'  >  /dev/null  2>&1 
r.null  $@_tl  setnull=0 
r.null  xlargeCityTimeRoad  setnull=0 
g. region  res=90 

r.cost  -r  start_rast=$@_tl  input=overlandTravelTime90  output=$@_t2  percent=100 
g. region  res=30 

r.mapcalc  xIargeCityTime=' if (isnull (xlargeCityTimeRoad) , $@_t2, xlargeCityTimeRoad) ' 
g. remove  $@_jtl,$@_t2 

#  Creates  a  slope  map  (in  degrees)  from  the  dem.  Assumes  the  dan  z  values  are  meters 
slope^att:  dem 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r. slope. aspect  dem  slope=slope_att_tmp 

r.mapcalc  'slope_att=round(slope_att_tmp) '  >  /dev/null  2>&1 
g. remove  slope_att_tmp 

#  Create  a  travel  cost  surface  for  travel  to  ramps  (road  cat  6) 
rampTimeRoad:  cross  intTravelTime30  othTravelTime30  otherroads  interstates 

@echo  1 ######################### 1 ;  echo  $@  ;  date 

#  Look  for  ramp  (cat  6)  cells  that  border  limited  access  roads  (cat  1) 
g. region  res=30 

r.mapcalc  $@__tl='if  (otherroads=6  &&  (\ 

interstates  [-1,-1]=!  ||  interstates  [-1,0]=1  ||  interstates  [-1, 1]  =1  ||  \ 
interstates  [  0,-1  ]=1  ||  interstates  [  0,1]=1  ||  \ 

interstates  [  1,-1  ]=1  ||  interstates  [  1,0]=1  ||  interstates  [  1, 1]  =1) ,  l,null  () )  '  > 
/dev/null  2>&1 

/data/FF/bin/r.cost.dks  in=othTravelTime30  m2=intTravelTime30  \ 
out=$@  xover=cross  start_rast=$@_tl  percent=100  \ 

>  /dev/null  2>&1 
g. remove  $@_tl 

#  Create  a  travel  cost  surface  for  travel  to  ramps 
ramp_att:  rampTimeRoad  overlandTravelTime90 

@echo  1 ######################### 1 ;  echo  $@  ;  date 
g. region  res=30 
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r.null  rarnpTimeRoad  null=0 
r.mapcalc  '$@_tl  =  max(  \ 
rairpTimeRoad  [-1 ,  -1] ,  \ 
rairpTimeRoad  [-1 , 0  ] ,  \ 
rairpTimeRoad  [-1 , 1  ] ,  \ 
rairpTimeRoad [  0,-1], \ 
rarrpTimeRoad  ,  \ 
rairpTimeRoad [  0,1],  \ 
rairpTimeRoad  [  1,-1],  \ 
rarrpTimeRoad  [  1,0],  \ 
rarrpTimeRoad  [  1,1])'  >  /dev/null  2>&1 
r.null  $@_tl  setnull=0 
r.null  rairpTimeRoad  setnull=0 
g. region  res=90 

r.cost  -r  start_rast=$@_tl  input=overlandTravelTime90  output=$@_t2  percent=100 
g. region  res=30 

r.mapcalc  rarnp_att=' if  (isnull  (rairpTimeRoad) ,  $@_t2, rairpTimeRoad)  ' 
g. remove  $@_tl,$@_t2 

#  Create  a  travel  cost  surface  for  travel  to  state  roads 
stateTimeRoad:  cross  intTravelTime30  othTravelTime30  otherroads 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

#  Look  for  state  roads 

r.mapcalc  $@__tl='if  (otherroads=3,  l,null  () )  '  >  /dev/null  2>&1 
/data/FF/bin/r.cost.dks  in=othTravelTime30  m2=intTravelTime30  \ 
out=$@  xover=cross  start_rast=$@_jtl  percent=100  \ 

>  /dev/null  2>&1 
g. remove  $@_tl 

#  Create  a  travel  cost  surface  for  travel  to  ramps 
staterd_att:  stateTimeRoad  overlandTravelTime90 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 
r.null  stateTimeRoad  null=0 
r.mapcalc  '$@_tl  =  max(  \ 
stateTimeRoad [-1,  -1] ,  \ 
stateTimeRoad [-1, 0] ,  \ 
stateTimeRoad [-1, 1] ,  \ 
stateTimeRoad [  0, -1 ] , \ 
stateTimeRoad  ,  \ 
stateTimeRoad [  0, 1 ] , \ 
stateTimeRoad [  1, -1] , \ 
stateTimeRoad [  1 , 0 ] , \ 
stateTimeRoad [  1,1]) '  >  /dev/null  2>&1 
r.null  $@_tl  setnull=0 
r.null  stateTimeRoad  setnull=0 
g. region  res=90 

r.cost  -r  start_rast=$@_tl  input=overlandTravelTime90  output=$@_t2  percent=100 
g. region  res=30 

r.mapcalc  staterd_att=' if (isnull (stateTimeRoad) , $@_t2, stateTimeRoad) ' 
g. remove  $@_tl,$@_t2 

#  Create  a  travel  cost  surface  for  travel  to  county  roads 
countyTimeRoad:  cross  intTravelTime30  othTravelTime30  otherroads 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

#  Look  for  county  roads 

r.mapcalc  $@_tl='if  (otherroads=4,  l,null  () )  '  >  /dev/null  2>&1 
/data/FF/bin/r.cost.dks  in=othTravelTime30  m2=intTravelTime30  \ 
out=$@_t2  xover=cross  start_rast=$@_tl  percent=100  \ 

>  /dev/null  2>&1 

r.mapcalc  '$@=if  ($@_t2=0,othTravelTime30, $@_t2)  ' 

#  g. remove  $@_tl,$@_t2 
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#  Create  a  travel  cost  surface  for  travel  to  ramps 

road_att:  countyTimeRoad  overlandTravelTime90 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 
r.null  countyTimeRoad  null=0 
r.mapcalc  '$@_tl  =  max(  \ 
countyTimeRoad [ -1 , -1 ] , \ 
countyTimeRoad [ -1 , 0 ] , \ 
countyTimeRoad [ -1 , 1 ] ,  \ 
countyTimeRoad [  0 , -1 ] , \ 
countyTimeRoad  ,  \ 
countyTimeRoad [  0 , 1 ] , \ 
countyTimeRoad [  1 , -1 ] , \ 
countyTimeRoad [  1 , 0 ] , \ 
countyTimeRoad [  1,1])'  >  /dev/null  2>&1 
r.null  $@_tl  setnull=0 
r.null  countyTimeRoad  setnull=0 
g. region  res=90 

r.cost  -r  start_rast=$@_tl  input=overlandTravelTime90  output=$@_t2  percent=100 
g. region  res=30 

r.mapcalc  road_att=' if (isnull (countyTimeRoad) , $@_t2, countyTimeRoad)  ' 
g. remove  $@_tl,$@_t2 

#  Create  a  travel  cost  surface  for  travel  to  intersections 

inter sectTimeRoad:  cross  intTravelTime30  othTravelTime30  otherroads  intersection 
@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

/data/FF/bin/r.cost.dks  in=othTravelTime30  m2=intTravelTime30  \ 
out=$@  xover=cross  start_rast=intersection  percent=100  \ 

>  /dev/null  2>&1 

#  Create  a  travel  cost  surface  for  travel  to  intersections 

intersect_att :  intersectTimeRoad  overlandTravelTime90 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

r.null  intersectTimeRoad  null=0 
r.mapcalc  '$@_tl  =  max(  \ 
intersectTimeRoad [-1, -1] ,  \ 
intersectTimeRoad  [-1 , 0] , \ 
intersectTimeRoad [-1 , 1 ] , \ 
intersectTimeRoad [  0 , -1 ] , \ 
intersectTimeRoad  , \ 
intersectTimeRoad  [  0 , 1 ] , \ 
intersectTimeRoad [  1, -1] , \ 
intersectTimeRoad  [  1 , 0 ] , \ 
intersectTimeRoad [  1,1])'  >  /dev/null  2>&1 
r.null  $@_tl  setnull=0 
r.null  intersectTimeRoad  setnull=0 
g. region  res=90 

r.cost  -r  start_rast=$@_tl  input=overlandTravelTime90  output=$@_t2  percent=100 
g. region  res=30 

r .mapcalc  intersect_att= ' if (isnull (intersectTimeRoad) , $@_t2 , intersectTimeRoad) ' 
g. remove  $@_tl,$@_t2 

#  Finds  state  (2)  and  county  (3)  road  intersections.  Create  a  map  of 

#  just  state  and  county  roads.  Expand  these  locations  by  one  cell  and 

#  then  thin  than  down  to  one  cell  (to  take  care  of  any  small  map  errors) . 

#  Then  find  cells  that  are  surrounded  by  more  than  two  cells  containing 

#  roads. 

intersection:  otherroads 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  '$@_tl=if  ( (otherroads=2)  ||  ( other r oads=3 ) ,  1,  0)'  >  /dev/null  2>&1 
r.thin  $@_tl  output=$@_t2 
r.null  $@  t2  null=0 
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r.mapcalc  ' intersection=if  ($@_t2=l, if (2<  \ 

($@_t2 [-1,-1]  +  $@_t2 [— 1, 0 ]  +  $@_t2 [-1,1]  +  \ 

$@_t2 [  0,-1]  +  $@_t2 [  0,1]  +  \ 

$@_t2 [  1,-1]  +  $@_t2 [  1,0]  +  $@_t2 [  1,1]), 1,0))'  >  /dev/ null  2>&1 
r.null  intersection  setnull=0 
g. remove  $@_tl,$@_t2 

#  Create  the  municipal  buffer  file  for  LEAM 
cityBuff:  rnunicipalboundary  boundary 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  'cityBuff=if (isnull (boundary) , 1, if (isnull (rnunicipalboundary) , 0, 1) ) '  >  /dev/null 

2>&1 

#  Create  a  surface  with  an  index  of  a  cell's  juxtaposition  with  residential  neighbors 
neighbor_att:  landcover 

@echo  ' ######################### ' ;  echo  $@  ;  date 

(  \ 

echo  TITLE  Find  dense  residential  areas  ;\ 

echo  MATRIX  29  ;\ 


echo 
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;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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echo 
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echo 
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echo 
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echo 
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echo 
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echo 
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1  ;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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;\ 

echo 
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;\ 
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;\ 
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;\ 
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;\ 
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;\ 
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;\ 
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;\ 

echo  DIVISOR  1  ;\ 
echo  TYPE  P  ;\ 

)  >  /tmp/filter$$ 

g. region  res=30 

#  select  urban  residential  areas 

r.mapcalc  $@_tl='if  (landcover=21, 1,  if  (landcover =2  2, 2, 0) )  ' 
r.mfilter  input=$@_tl  output=neighbor_att  filter=/trnp/filter$$ 
g. remove  $@_tl 
rm  /trrp/filter$$ 

#  Create  a  distance  surface  to  the  nearest  cell  containing  water. 
water_att:  landcover 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 
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r.mapcalc  $@_tl='if  (landcover=ll  ||  landcover=91  ||  landcover=92, 1,0)  '  >  /dev/null 

2>&1 

r.null  $@_tl  setnull=0 

r. buffer  $@_tl  output=$@_t2  distances=30, 60, 90, 120 

r.mapcalc  water_att='if  (isnull  ($@__t2) ,  127,  ($@_t2-l)  *30)  '  >  /dev/null  2>&1 
g. remove  $@_tl,$@_t2 

#  Create  a  distance  surface  to  the  nearest  cell  containing  forest  - 
forest_att:  landcover 

@echo  1 ######################### 1 ;  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  $@_tl='if  (landcover=41  ||  landcover=42  ||  \ 
landcover=43  ||  landcover=91, 1, 0)  '  >  /dev/null  2>&1 
r.null  $@_tl  setnull=0 

r. buffer  $@_tl  output=$@_t2  distances=30, 60, 90, 120 

r.mapcalc  forest_att=' if (isnull ($@_t2) , 127, ($@_t2-l)  *30) '  >  /dev/null  2>&1 
g. remove  $@_tl,$@_t2 

###################################################################### 

###################################################################### 

#  Travel  time  maps  generated  below  at  the  30  and  90  meter  resolution 
###################################################################### 
###################################################################### 

#  Creates  a  travel  time  across  90  meter  cells  by  choosing  the  shortest  time 

#  from  among  the  9  associated  30  meter  cells.  It  preserves  travel  time  along 

#  road  cat  1  (limited  access)  and  preserves  the  highway  buffer 
overlandTravelTime90 :  overlandTravelSpeed90 

@echo  1 ######################### 1 ;  echo  $@  ;  date 
g. region  res=90 

#  90m/cell  *  lmile/1609meter  *  60min/hr  =  3. 35538min. mile/cell. hr 

#  3. 35538min. mile/cell. hr  /  x  miles/hr  =  Y  min/cell 

r .mapcalc  ' overlandTravelTime90=if (overlandTravelSpeed90>0, 3 . 35538/overlandTravelSpeed90, 
nullQ)  '  >  /dev/null  2>&1 

#  Computes  the  time  required  to  traverse  a  cell  in  the  N-S  or  E-W  axes 
roadTime30:  roadSpeed30 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

#  30m/cell  *  lmile/1609meter  *  60min/hr  =  1 .11846min. mile/cell. hr 

#  1.11846min. mile/cell. hr  /  x  miles/hr  =  Y  min/cell 

r.mapcalc  ' roadTime30=if (roadSpeed30>0, 1 .11846/roadSpeed30,  null())'  >  /dev/null  2>&1 
#travelTime30 :  travelSpeed30 

#  @echo  1 ######################### 1 ;  echo  $@  ;  date 

#  g. region  res=30 

#  #  30m/cell  *  lmile/1609meter  *  60min/hr  =  1 .11846min. mile/cell. hr 

#  #  1 . 11846min. mile/cell. hr  /  x  miles/hr  =  Y  min/cell 

#  r.mapcalc  'travelTime30=if (travelSpeed30>0, 1 .11846/travelSpeed30,  null())'  >  /dev/null 
2>&1 

othTravelTime30 :  othTravelSpeed30 

@echo  1 ######################### 1 ;  echo  $@  ;  date 
g. region  res=30 

#  30m/cell  *  lmile/1609meter  *  60min/hr  =  1 .11846min. mile/cell. hr 

#  1.11846min. mile/cell. hr  /  x  miles/hr  =  Y  min/cell 

r.mapcalc  'othTravelTime30=if (othTravelSpeed30>0, 1 .11846/othTravelSpeed30,  null())'  > 
/dev/null  2>&1 

intTravelTime30 :  intTravelSpeed30 

@echo  1 ######################### 1 ;  echo  $@  ;  date 
g. region  res=30 

#  30m/cell  *  lmile/1609meter  *  60min/hr  =  1 .11846min. mile/cell. hr 

#  1.11846min. mile/cell. hr  /  x  miles/hr  =  Y  min/cell 

r.mapcalc  'intTravelTime30=if (intTravelSpeed30>0, 1 .11846/intTravelSpeed30,  null())'  > 
/dev/null  2>&1 
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overlandTravelSpeed90 :  hwyBuff  landTravelSpeed30 

@echo  1 ######################### 1 ;  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  overlandTravelSpeed90='eval (  \ 
highland  =  max(  \ 

landTravelSpeed30 [-1, -1] , \ 
landTravelSpeed30 [-1 , 0] , \ 
landTravelSpeed30 [-1, 1] , \ 
landTravelSpeed30 [  0 , -1 ] , \ 
landTravelSpeed30  , \ 
landTravelSpeed30 [  0 , 1 ] , \ 
landTravelSpeed30 [  1, -1] , \ 
landTravelSpeed30 [  1 , 0] , \ 
landTravelSpeed30 [  1,1]),  \ 
haveBuffer  =  if (  \ 

hwyBuff  [-1,-1]  =1  ||  hwyBuff  [-1,0]  =1  ||  hwyBuff  [-1,1]  =1  ||  \ 

hwyBuff  [  0,-l]=l  ||  hwyBuff=l  I  I  hwyBuff  [  0,1]=1  ||  \ 

hwyBuff [  1,-1]=1  ||  hwyBuff [  1,0]=1  I  I  hwyBuff[  1,1]  =1,  1,  0),  \ 

if (haveBuffer, 0, highland) ) '  >  /dev/null  2>&1 

intTravelSpeed30 :  interstates 

@echo  1 ######################### 1 ;  echo  $@  ;  date 

g. region  res=30 

r.mapcalc  intTravelSpeed30='\ 

if  (interstates=l,70,0)  '  >  /dev/null  2>&1 

othTravelSpeed30 :  otherroads 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 
r.mapcalc  othTravelSpeed30='\ 
if  (otherroads=l,  70,  \ 
if  (otherroads=2, 50,  \ 
if  (otherroads=3, 40,  \ 
if  (otherroads=4, 30,  \ 
if  (otherroads=5, 30,  \ 
if  (otherroads=6, 30,  \ 
if  (otherroads=7, 10,  \ 

if  (otherroads=8, 20, 0) )))))))  '  >  /dev/null  2>&1 

landTravelSpeed30:  landcover  noGrowth 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 
r .mapcalc  landTravelSpeed30= ' \ 
if  (noGrowth=l,  0,  \ 
if  (landcover=21  |  |  \ 
landcover=22  |  |  \ 
landcover=23  ,.5,  \ 
if  (landcover=ll  |  |  \ 
landcover=12  |  |  \ 
landcover=91  |  |  \ 

landcover=92,  0,  1) ) )  '  >  /dev/null  2>&1 

#  r.mapcalc  landTravelSpeed30= ' \ 

#  if  (landcover=21, 10,  \ 

#  if  (landcover=22, 10,  \ 

#  if  (landcover=23, 10,  null  () ) ) )  '  \  >  /dev/null  2>&1 

# 

#travelablelandSpeed:  landcover 

#  r.mapcalc  travelablelandSpeed='  \ 

#  if (noGrowth ! =1  &&  \ 

#  landcover! =11  &&  \ 

#  landcover ! =12  &&  \ 

#  landcover! =91  &&  \ 

#  landcover! =92,  1,  0) '  >  /dev/null  2>&1 
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#  Creates  a  map  of  1-cell  wide  buffers  along  limited  access  highways 

#  except  where  other  roads  intersect.  Problems: 

#  -  No  distinction  between  bridges  and  true  intersections 

#  -  Parallel  roads  in  the  next  cell  are  considered  access  points 

hwyBuff:  interstates  otherroads 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=30 
r.mapcalc  hwyBuff='\ 

if (interstates  =  1  |  f  otherroads  =6,  0, \ 

if  (interstates  [-1,-1]  =1  ||  interstates  [-1, 0]  =1  ||  interstates  [-1,1]=1  ||  \ 
interstates  [  0,-1  ]=1  |  |  interstates  [  0,1]=1  |  |  \ 

interstates[  1,-1]=1  ||  interstates[  1,0]  =1  II  interstates[  1,1]=1,  1,  0))'  > 
/dev/null  2>&1 

###################################################################### 

###################################################################### 

#  The  following  instructions  generate  index  (0-1)  maps  for  each 

#  attractiveness  factor  using  the  following  approach: 

#  1  -  cross-list  the  attractiveness  level  to  starting  residential  patterns 

#  and  potential  developable  areas. 

#  2  -  associate  attractiveness  level  with  percent  developed 

#  3  -  convert  attractiveness  map  to  percent  developed  (0-1  index) 

# 

#  This  is  accomplished  with  these  steps: 

#  1  -  Find  high  value  in  attractiveness  map  (r. describe  -r) 

#  2  -  Divide  the  map  by  this  value;  multiply  by  10  and  round  to  create 

#  10  categories  (r.mapcalc) 

#  3  -  Cross-reference  the  categories  with  the  number  of  residential  and 

#  developable  cells  (r.coin) 

#  4  -  Process  the  results  to  create  a  "graph"  (awk) 

#  5  -  Apply  the  graph  to  the  attractiveness  map  to  create  the  index  map 
###################################################################### 
###################################################################### 

cities_att2 :  cities_att 

r.mapcalc  cities_att2='log(if  (row()=l  &&  col  0=1, 0,cities_att) )  ' 

developable:  landcover  noGrowth  hwyBuff 

@echo  '#########################';  echo  $@  ;  date 
r.mapcalc  developable=  \ 

'  if  (hwyBuff  |  |  noGrowth  |  |  landcover=ll  |  |  landcover=24  |  |  \ 
landcover=91  |  |  landcover=92  |  |  landcover=85,  0,  \ 
if  (landcover=21  ||  landcover=22  ||  landcover=24,  2,  1))'  \ 

>  /dev/null  2>&1 

. ,/LEAM/XXcitiesGraph:  . . /LEAM/awkscript  developable  cities_att2 
@echo  '#########################';  echo  $@  ;  date 
@  rm  -f  . . /LEAM/XXcitiesHigh 
@  mkdir  -p  . . /LEAM 
##  Normalize  the  map 

#  find  highest  value  in  map 

r. describe  cities_att2  -r  -i  -q  |  grep  thru  |  sed  -e  's/.*  //'  >  . . /IEAM/XXcitiesHigh 

#  divide  map  by  highest  value  -  multiply  by  some  number  (e.g.,  20) 

r.mapcalc  XXcitiesGraph_tl=' round ( (20  *  cities_att2) / ' 'cat  ../LEAM/XXcitiesHigh'')'  > 
/dev/null  2>&1 

#  generate  awkscript  for  processing  r. stats  output 

sed  -e  s:MULT:'cat  ../LEAM/XXcitiesHigh':  ../LEAM/awkscript  >  . . /LEAM/awk-cities 

#  apply  graph  to  attractiveness  map  to  create  index  map 

r. stats  -c  XXcitiesGraph_tl, developable  |  awk  -f  ../LEAM/awk-cities  > 

. . /LEAM/XXcitiesGraph 

g. remove  XXcitiesGraph_tl 

cities_prob:  cities_att2  ../LEAM/XXcitiesGraph 
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r.mapcalc  ' cities_prob^graph (float (cities_att2) ,  ''cat  . ,/LEAM/XXcitiesGraph' ' ) '  > 

/dev/ null  2>&1 

(echo  0.0  white;  echo  .01  255  240  150;  echo  .05  yellow;  echo  .4  green;  \ 
echo  1.0  0  100  0;echo  end)  |  r. colors  $@  color=rules  -q 

. ,/LEAM/XXroadGraph:  ../LEAM/awkscript  developable  road_att2 
@echo  ' ######################### ' ;  echo  $@  ;  date 
@  rm  -f  . . /LEAM/XXroadHigh 
@  mkdir  -p  . . /LEAM 
##  Normalize  the  map 

#  find  highest  value  in  map 

r. describe  road_att2  -r  -i  -q  |  grep  thru  |  sed  -e  's/.*  //'  >  ../LEAM/XXroadHigh 

#  divide  map  by  highest  value  -  multiply  by  some  number  (e.g.,  20) 

r.mapcalc  XXroadGraph_tl= ' round ( (20  *  road_att2) / ' 'cat  . ./LEAM/XXroadHigh' 1 ) '  >  /dev/null 

2>&1 

#  generate  awkscript  for  processing  r. stats  output 

sed  -e  s:MULT:'cat  . ./LEAM/XXroadHigh' :  ../LEAM/awkscript  >  . . /LEAM/awk-road 

#  apply  graph  to  attractiveness  map  to  create  index  map 

r. stats  -c  XXroadGraph_tl, developable  |  awk  -f  . . /LEAM/awk-road  >  . . /LEAM/XXroadGraph 
g. remove  XXroadGraph_tl 

road_att2:  road_att 

r.mapcalc  ’road_att2=log(10*road_att)  1 

road_prob:  road_att2  developable  ../LEAM/XXroadGraph 

r.mapcalc  1 road_prob=graph (float (road_att2) ,  ’'cat  . ./LEAM/XXroadGraph' ') '  >  /dev/null 

2>&1 

(echo  0.0  white;  echo  .01  255  240  150;  echo  .05  yellow;  echo  .4  green;  \ 
echo  1.0  0  100  0;echo  end)  |  r. colors  $@  color=rules  -q 

.  . /LEAM/XXrampGraph:  . .  /LEAM/awkscript  developable  rarrp_att2 
@echo  ' ######################### ' ;  echo  $@  ;  date 
@  rm  -f  . .  /LEAM/XXrampHigh 
@  mkdir  -p  . . /LEAM 
##  Normalize  the  map 

#  find  highest  value  in  map 

r. describe  ramp_att2  -r  -i  -q  |  grep  thru  |  sed  -e  's/.*  //'  >  . . /LEAM/XXrampHigh 

#  divide  map  by  highest  value  -  multiply  by  some  number  (e.g.,  20) 

r.mapcalc  XXrampGraph_tl= ' round ( (20  *  ramp_att2) / ' 'cat  ../LEAM/XXrampHigh'')'  >  /dev/null 

2>&1 

#  generate  awkscript  for  processing  r. stats  output 

sed  -e  s:MULT:'cat  ../LEAM/XXrampHigh':  . .  /LEAM/awkscript  >  . .  /LEAM/awk-ramp 

#  apply  graph  to  attractiveness  map  to  create  index  map 

r.  stats  -c  XXrampGraph_tl,  developable  |  awk  -f  ../LEAM/awk-ramp  >  . . /LEAM/XXrairpGraph 
g. remove  XXrarrpGraph_tl 

ramp_att2:  rairp^att 

r.mapcalc  'ramp_att2=log(10*rarnp_att)  ' 

ramp_prob:  ramp_att2  developable  . . /LEAM/XXrampGraph 

r.mapcalc  '  ramp_prob=graph  (float  (ramp_att2) ,  ''cat  ../LEAM/XXrairpGraph'')'  >  /dev/null 

2>&1 

(echo  0.0  white;  echo  .01  255  240  150;  echo  .05  yellow;  echo  .4  green;  \ 
echo  1.0  0  100  0;echo  end)  |  r. colors  $@  color=rules  -q 

. ,/LEAM/XXintersectGraph:  ../LEAM/awkscript  developable  inter sect_att2 
@echo  ' ######################### ' ;  echo  $@  ;  date 
@  rm  -f  . . /LEAM/XXintersectHigh 
@  mkdir  -p  . . /LEAM 
##  Normalize  the  map 

#  find  highest  value  in  map 

r. describe  intersect_att2  -r  -i  -q  |  grep  thru  |  sed  -e  's/.*  //'  > 

. . /LEAM/XXintersectHigh 

#  divide  map  by  highest  value  -  multiply  by  some  number  (e.g.,  20) 
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r.mapcalc  XXintersectGraph_tl= ' round ( (20  *  inter sect_att2) /' 'cat 
. ,/LEAM/XXintersectHigh' ' ) '  >  /dev/null  2>&1 

#  generate  awkscript  for  processing  r. stats  output 

sed  -e  s:MULT:'cat  . ,/LEAM/XXintersectHigh' :  . . /LEAM/awkscript  >  . . /LEAM/awk- intersect 

#  apply  graph  to  attractiveness  map  to  create  index  map 

r. stats  -c  XXintersectGraph_tl, developable  |  awk  -f  . . /LEAM/awk-intersect  > 

. . /LEAM/XXintersectGraph 

g. remove  XXintersectGraph_tl 

intersect_att2 :  intersect_att 

r.mapcalc  intersect_att2='log(10*intersect_att) '  ; 

inter sect_prob:  intersect_att2  developable  ../LEAM/XXintersectGraph 

r.mapcalc  ' intersect_prob=graph (float (intersect_att2) ,  ''cat  . ./LEAM/XXintersectGraph' ' ) ' 
>  /dev/null  2>&1 

(echo  0.0  white;  echo  .01  255  240  150;  echo  .05  yellow;  echo  .4  green;  \ 
echo  1.0  0  100  0;echo  end)  |  r. colors  $@  color=rules  -q 

. ,/LEAM/XXforestGraph:  . . /LEAM/ awkscript  developable  forest_att 

@echo  ' ######################### ' ;  echo  XXforestGraph  ;  date 
@  rm  -f  . . /LEAM/XXforestHigh 
@  mkdir  -p  . .  /LEAM 
##  Normalize  the  map 

#  find  highest  value  in  map 

r. describe  forest_att  -r  -i  -q  |  grep  thru  |  sed  -e  's/.*  //'  >  ../LEAM/XXforestHigh 

#  divide  map  by  highest  value  -  multiply  by  some  number  (e.g.,  20) 

r.mapcalc  XXforestGraph_tl=' round ( (20  *  forest_att) / ' 'cat  . ./LEAM/XXforestHigh' ’ ) ’  > 
/dev/null  2>&1 

#  generate  awkscript  for  processing  r. stats  output 

sed  -e  s:MULT:'cat  ../LEAM/XXforestHigh':  . . /LEAM/awkscript  >  . . /LEAM/awk- forest 

#  apply  graph  to  attractiveness  map  to  create  index  map 

r. stats  -c  XXforestGraph_tl, developable  |  awk  -f  . . /LEAM/awk- f orest  > 

. . /LEAM/XXforestGraph 

g. remove  XXforestGraph_tl 

forest_prob:  forest_att  developable  ../LEAM/XXforestGraph 

r.mapcalc  ' forest_prob=graph (float (forest_att) ,  ''cat  . . /IEftM/XXforestGraph' 1 ) ’  > 
/dev/null  2>&1 

(echo  0.0  white;  echo  .01  255  240  150;  echo  .05  yellow;  echo  .4  green;  \ 
echo  1.0  0  100  0;echo  end)  |  r. colors  $@  color=rules  -q 

. ,/LEflM/XXwaterGraph:  . . /LEAM/awkscript  developable  water_att 

@echo  ' ######################### ' ;  echo  XXwaterGraph  ;  date 
@  rm  -f  . . /LEflM/XXwaterHigh 
@  mkdir  -p  . . /IEAM 
##  Normalize  the  map 

#  find  highest  value  in  map 

r. describe  water_att  -r  -i  -q  |  grep  thru  I  sed  -e  's/.*  //'  >  . . /LEAM/XXwaterHigh 

#  divide  map  by  highest  value  -  multiply  by  some  number  (e.g.,  20) 

r.mapcalc  XXwaterGraph_tl=' round ( (20  *  water_att) /' 'cat  ../LEAM/XXwaterHigh'’)'  > 
/dev/null  2>&1 

#  generate  awkscript  for  processing  r. stats  output 

sed  -e  s:MULT:'cat  ../LEAM/XXwaterHigh':  ../LEAM/awkscript  >  . . /LEAM/awk-water 

#  apply  graph  to  attractiveness  map  to  create  index  map 

r. stats  -c  XXwaterGraph_tl, developable  |  awk  -f  ../LEAM/awk-water  >  . . /LEAM/XXwaterGraph 
g. remove  XXwaterGraph_tl 

water_prob:  water_att  developable  ../LEAM/XXwaterGraph 

r.mapcalc  ’ water_prob=graph ( float (water_att ) ,  ’'cat  ../LEAM/XXwaterGraph'1)'  >  /dev/null 

2>&1 

(echo  0.0  white;  echo  .01  255  240  150;  echo  .05  yellow;  echo  .4  green;  \ 
echo  1.0  0  100  0;echo  end)  |  r. colors  $@  color=rules  -q 

. ,/LEAM/XXneighborGraph:  ../LEAM/awkscript  developable  neighbor_att 

@echo  ' ######################### ' ;  echo  XXneighborGraph  ;  date 
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@  rm  -f  . . /LEAM/XXneighborHigh 
@  mkdir  -p  . . /LEAM 
##  Normalize  the  map 

#  find  highest  value  in  map 

r. describe  neighbor_att  -r  -i  -q  |  grep  thru  |  sed  -e  's/.*  //'  >  ../LEAM/XXneighborHigh 

#  divide  map  by  highest  value  -  multiply  by  some  number  (e.g.,  20) 

r.mapcalc  XXneighborGraph_tl=' round ( (20  *  neighbor_att) / ' 'cat  . ./LEAM/XXneighborHigh' ' ) ' 

>  /dev/ null  2>&1 

#  generate  awkscript  for  processing  r. stats  output 

sed  -e  s:MULT:'cat  ../LEAM/XXneighborHigh':  .. /LEAM/ awkscript  >  . . /LEAM/awk-neighbor 

#  apply  graph  to  attractiveness  map  to  create  index  map 

r. stats  -c  XXneighborGraph_tl, developable  |  awk  -f  ../LEAM/awk-neighbor  > 

. . /LEAM/XXneighborGraph 

g. remove  XXneighborGraph_tl 

neighbor_prob :  neighbor_att  developable  . . /LEAM/XXneighborGraph 

r.mapcalc  ' neighbor_prob=graph( float ( neighbor_att ) ,  ''cat  ../LEAM/XXneighborGraph'')'  > 
/dev/null  2>&1 

(echo  0.0  white;  echo  .01  255  240  150;  echo  .05  yellow;  echo  .4  green;  \ 
echo  1.0  0  100  0;echo  end)  |  r. colors  $@  color=rules  -q 

. ,/LEAM/XXslopeGraph:  .. /LEAM/ awkscript  developable  slope_att 
@echo  ' ######################### ' ;  echo  $@  ;  date~ 

@  rm  -f  . .  /LEAM/XXslopeHigh 
@  mkdir  -p  . . /LEAM 
##  Normalize  the  map 

#  find  highest  value  in  map 

r. describe  slope_att  -r  -i  -q  |  grep  thru  |  sed  -e  's/.*  //'  >  ../LEAM/XXslopeHigh 

#  divide  map  by  highest  value  -  multiply  by  some  number  (e.g.,  20) 

r.mapcalc  XXslopeGraph_tl=' round ( (20  *  slope_att) / ' ' cat  ../LEAM/XXslopeHigh'1)'  > 
/dev/null  2>&1 

#  generate  awkscript  for  processing  r. stats  output 

sed  -e  s:MULT:'cat  ../LEAM/XXslopeHigh':  .. /LEAM/ awkscript  >  . . /LEAM/awk-slope 

#  apply  graph  to  attractiveness  map  to  create  index  map 

r. stats  -c  XXslopeGraph_tl, developable  |  awk  -f  ../LEAM/awk-slope  >  . . /LEAM/XXslopeGraph 
g. remove  XXslopeGraph_tl 

slope_prob:  slope_att  developable  ../LEAM/XXslopeGraph 

r.mapcalc  ’ slope_prob=graph ( float ( slope_att ) ,  ’'cat  ../LEAM/XXslopeGraph'1)'  >  /dev/null 

2>&1 

(echo  0.0  white;  echo  .01  255  240  150;  echo  .05  yellow;  echo  .4  green;  \ 
echo  1.0  0  100  0;echo  end)  |  r. colors  $@  color=rules  -q 

. ,/LEAM/XXstaterdGraph:  . . /LEAM/awkscript  developable  staterd_att2 
@echo  ' ######################### ' ;  echo  $@  ;  date 
@  rm  -f  . . /LEAM/XXstaterdHigh 
@  mkdir  -p  . . /LEAM 
##  Normalize  the  map 

#  find  highest  value  in  map 

r. describe  staterd_att2  -r  -i  -q  |  grep  thru  |  sed  -e  's/.*  //'  >  ../LEAM/XXstaterdHigh 

#  divide  map  by  highest  value  -  multiply  by  some  number  (e.g.,  20) 

r.mapcalc  XXstateRoadGraph_tl=' round ( (20  *  staterd_att2) /' 'cat  ../LEAM/XXstaterdHigh'')' 

>  /dev/null  2>&1 

#  generate  awkscript  for  processing  r. stats  output 

sed  -e  s:MULT:'cat  ../LEAM/XXstaterdHigh':  ../LEAM/awkscript  >  . . /LEAM/awk-staterd 

#  apply  graph  to  attractiveness  map  to  create  index  map 

r. stats  -c  XXstateRoadGraph_tl, developable  |  awk  -f  ../LEAM/awk-staterd  > 

. . /LEAM/XXstaterdGraph 

g. remove  XXstatePoadGraph_tl 

staterd_att2 :  staterd_att 

r.mapcalc  staterd_att2='log(10*staterd_att) ' 

staterd_prob:  staterd_att2  ../LEAM/XXstaterdGraph 
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r.mapcalc  ' staterd_prob=graph (float (staterd_att2) ,  ''cat  . . /LEAM/XXstaterdGraph' ' ) '  > 
/dev/null  2>&1 

(echo  0.0  white;  echo  .01  255  240  150;  echo  .05  yellow;  echo  .4  green;  \ 
echo  1.0  0  100  0;echo  end)  |  r. colors  $@  color=rules  -q 

###################################################################### 

#  This  combines  residential  attractors 

#  into  an  overall  residential  attractor  map 

#  residential_att_mult :  developable  road_prob  slope_prob  ramp_prob  \ 

#  staterd_prob  forest_prob  water_prob  neighbor_prob\ 

#  cities_prob  intersect_prob  noGrowth  landcover  boundary 

#  @echo  T#########################';  echo  $@  ;  date 

#  r.mapcalc  residential_att_mult=eval ' (  \ 

#  developable=if (developable !=0, 1,0),  \ 

#  developable  *  road_prob  *  slope_prob  *  intersect_prob  *  \ 

#  rarrp_prob  *  staterd_prob  *  forest_prob  *  water_prob  *  \ 

#  neighbor_prob  *  cities_prob  ) '  >  /dev/null  2>&1 

residential_att:  developable  road_prob  slope_prob  ramp_prob  \ 
staterd_prob  forest_prob  water_prob  neighbor_prob\ 
cities_prob  intersect_prob  noGrowth  landcover  boundary 
@echo  T#########################';  echo  $@  ;  date 
r.mapcalc  residential_att=eval ' (  \ 
developable=if (developable>0,l,0),  \ 
developable  *  (road_prob  *  ${road_prob_MULT}  +  \ 
slope_prob  *  ${ slope  jorob_MULT}  +  \ 
inter sect_prob  *  ${intersect_prob_MULT}  +  \ 
ramp_prob  *  ${rarnp_prob_MULT}  +  \ 
staterd_prob  *  ${staterd_prob_MULT}  +  \ 
forest_prob  *  ${forest_probJdLT}  +  \ 
water_prob  *  ${water_prob_MQLT}  +  \ 
cities_prob  *  ${cities_prob_MULT})  /  \ 

(${intersect_prob_MJLT}  +  ${road_prob_MULT}  +  \ 

${slope_prob_MJLT}  +  \ 

${rarnp_prob_MULT}  +  ${staterd_prob_MULT}  +  \ 

${forest_prob_MULT}  +  ${water_prob_MQLT}  +  \ 

${cities_prob_MULT} ) ) '  >  /dev/null  2>&1 

. ,/LEflM/XXresidentialGraph:  . . /LEAM/awkscript  developable  residential_att 
@echo  ' ######################### ' ;  echo  $@  ;  date 
@  rm  -f  . . /LEflM/XXattract_res_High 
@  mkdir  -p  . . /LEAM 
##  Normalize  the  map 

#  find  highest  value  in  map 

r. describe  residential_att  -r  -q  |  sed  -e  's/.*-//'  -e  's/  *$$//'  > 

. . /LEAM/XXattract_res_High 

#  divide  map  by  highest  value  -  multiply  by  some  number  (e.g.,  20) 
r.mapcalc  ResidentialGraph_tl= ' round ( (20  *  residential_att) / ' ' cat 

. ,/LEAM/XXattract_res_High' ' ) '  >  /dev/null  2>&1 

#  generate  awkscript  for  processing  r. stats  output 

sed  -e  s:MULT:'cat  . ,/IEAM/XXattract_res_High' :  ../LEAM/awkscript  >  ../LEAM/awk- 
residential 

#  apply  graph  to  attractiveness  map  to  create  index  map 

r. stats  -c  ResidentialGraph_tl, developable  |  awk  -f  . . /LEAM/awk-residential  > 

. . /LEAM/XXresidentialGraph 

g. remove  ResidentialGraph_tl 

residential_prob:  residential_att  . . /LEAM/XXresidentialGraph 

r.mapcalc  ’residential_prob=graph(float(residential__att) ,  ’'cat 
. . /LEAM/XXresidentialGraph' ') '  >  /dev/null  2>&1 

(echo  0.0  white;  echo  .01  255  240  150;  echo  .05  yellow;  echo  .4  green;  \ 
echo  1.0  0  100  0;echo  end)  |  r. colors  $@  color=rules  -q 

#  Evolve  Landscape 

futureLanduse :  residential_prob  landcover 
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r.mapcalc  newLanduse=if ' ( (residential_prob  >  rand (0. 0, ' 'r. describe  -r  residential_prob  | 
sed  -e  's/. *-//'''))  &&  (50>rand (0,100) ) ,21, null () ) '  >  /dev/null  2>&1 

r.mapcalc  futureLanduse='if  (landcover=21  ||  landcover=22, 1, if  (newLanduse=21,2, 0) )  '  > 
/dev/null  2>&1 

(echo  0  white;  echo  1  grey;  echo  2  red;  echo  end)  |  \ 
r. colors  $@  color=rules  -q 


noiseTolerance:  futurelanduse 

r.mapcalc  $@_tl='if (isnull (futureLanduse) ,0,1) '  >  /dev/null  2>&1 
r.mapcalc  $@_t2='  \ 

$@_tl [-2, -2] +$@_tl [-2, -1] +$@_tl [-2,0] +$@_tl [-2,1] +$@_tl [-2,2]+  \ 

$@_tl [-1, -2]+$@_tl [-1, -1] +$@_tl [-1,0] +$@_tl [-1,1] +$@_tl [-1,2]+  \ 

$@_tl[  0,-2]  +$@_tl  [  0,  —1]  +$@_tl;[  0,0]+$@_tl[  0,l]+$@_tl[  0,2]+  \ 

$@_tl[  1,-2]  +$@__tl [  l,-l]+$@_tl[  l,0]+$@_tl[  l,l]+$@_tl[  1,2]+  \ 

$@_tl[  2,-2] +$@_tl [  2,-1] +$@_tl [  2,0] +$@_tl [  2,1] +$@_tl [  2,2]'  >  /dev/null  2>&1 

g. region  res=150 

r. decay  $@_t2  out=noiseTolerance  dist=1000  ccmp=.01  per=100 
g. region  res=30 


altf utureLanduse :  residential_att 

#  x=$ (shell  r. stats  -c  residential_att  |  sed  -e  's/-[A  ]*/*/'  -e  's/$$/+\\/'  -e 
's/A\*.*/0/'  |  be  ) ;  \ 

#  echo  $$x  ;  \ 

r.mapcalc  futurelanduse=if ' (residential_att  *  ATTRACTOR_RES  >\ 
rand  (0 . 0, 1 . 0) , 21, null ( ) ) '  >  /dev/null  2>&1 


#  Removes  all  the  temporary  files  used  in  the  creation  of  the  LEAM  inputs 
clean: 

g. remove  rast=${Otherfiles] 


#  Removes  all  the  temporary  and  LEAM  input  files 
veryclean:  clean 

g. remove  rast=${LEfiles] 
g. remove  rast=$[Indexfiles] 
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Appendix  B:  Acquiring  and  Processing 
Required  GIS  Data 

The  following  provides  a  step-by-step  cookbook  for  acquiring  and  process¬ 
ing  required  GIS  data.  Remember,  however,  that  it  is  the  end  result  that  is 
important,  not  the  steps  to  get  there.  These  steps  are  provided  as  an  ap¬ 
proached  that  has  worked,  but— due  to  changes  in  available  data,  differ¬ 
ences  among  GIS  software,  your  hardware,  and  available  skills— different 
and  perhaps  better  approaches  will  be  used. 

The  steps  described  below  are  partly  accomplished  with  ESRI  GIS  soft¬ 
ware  and  partly  with  the  public  domain  GRASS  software. 

ESRI  Software  Steps 

Acquire  land  use  data  for  the  area  of  interest 

Acquire  GIS  data  layers  from  available  Internet  available  sources. 

•  Land  Cover:  USGS  NLCD  maps 

•  Digital  Elevation:  USGS  DEM 

•  Road  Network:  Census  Bureau  Tiger  files 

•  Property  Ownership:  USGS,  DoD,  BLM 

•  Floodplain  Areas:  FEMA 

•  County. 

Re-project  all  Data  into  the  Same  Equal  Area  Projection 

While  GRASS  can  process  vector  and  raster  GIS  data  in  virtually  any  pro¬ 
jection  and  coordinate  system,  it  cannot  simultaneously  work  with  maps 
that  are  in  different  projections.  Therefore  ArcGIS  is  used  to  re-project  all 
data  into  an  equal-area  coordinate  system  -  preferably  Universal  Trans¬ 
verse  Mercator  (UTM).  Also,  choose  a  minimal  bounding  box  to  be  your 
study  area  into  which  the  maps  will  be  re-projected. 

Process  the  Road  Data 

Because  GRASS  generally  works  with  single  category  values  (rather  than  a 
table  of  attributes),  vector  data  imported  into  GRASS  must  be  associated 
with  useful  category  values.  The  CFCC2  category  codes  associated  with 
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roads  are  string  values  such  as  “Ai”,  “A2”,  “A3”,  etc.  These  must  be  trans¬ 
formed  to  integer  values  by  dropping  the  “A.”  To  accomplish  this,  edit  the 
attributes  table  and  add  a  new  field  that  you  will  call  CFCC3  as  a  short  in¬ 
teger.  In  the  attribute  table  editor  for  CFCC3,  click  on  the  String  radio  but¬ 
ton.  Then  enter  the  text  in  the  box:  “Right  (CFCC2, 1)”  and  run  this  re¬ 
quest.  That  should  create  the  road  classes. 

Create  the  “No  Growth”  map. 

Use  the  Federal  lands  in  the  study  area  to  generate  the  no  growth  Shape- 
file.  Values  of  o  represent  areas  where  growth  is  allowed;  a  value  of  1  is 
used  for  areas  where  growth  is  restricted.  These  lands  are  designated  no 
growth  because  the  government  owns  them  and  residential  growth  will  not 
occur  there. 

Prepare  for  Copying  to  Your  GRASS  Environment 

At  this  point,  you  should  have  the  following  maps  all  in  the  same  UTM  co¬ 
ordinate  system: 


Code 

Name 

Description 

roads.shp 

Roads 

Vector  shape  file 

nogrowth. shp 

No  Growth 

Vector  shape  file 

county.shp 

Counties  in  study 

Vector  shape  file 

mnbnd 

Municipal  areas 

City  boundary  shape  file 

landcover 

Landcover 

Raster  grids  with  NLCD  categories  0-99 

DEM 

DEM 

Raster  grids  with  elevations  in  meters 

GRASS  Software  Steps 

The  Shapefiles  and  Rasters  must  be  copied  to  the  LINUX  environment  for 
processing  by  GRASS.  For  those  using  GRASS  under  Linux  this  will  re¬ 
quire  copying  the  maps  from  a  Windows  machine  to  the  LINUX  machine. 
(For  those  running  GRASS  under  CYGWIN  that  is  operating  under  Win¬ 
dows,  you  may  be  able  to  access  the  same  files  directly  without  copying.) 
Visit  your  local  system  administrator  for  help  using  your  systems,  net¬ 
work,  and  software. 

ESRI  vector  Shapefiles  should  be  copied  over  to  UNIX  as  a  set  of  six  files. 
These  files  will  have  the  same  name  but  different  extensions  (.dbf,  .prj, 
.sbn,  .sbx,  .shp,  .shx). 


Raster  grids  should  be  copied  over  to  UNIX  as  folders. 
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Create  New  GRASS  Location 

Start  GRASS  and  choose  to  create  a  new  location  if  this  is  the  first  time  for 
this  particular  study  area.  For  a  location,  choose  a  name  without  spaces 
(e.g.,  FortBragg),  and  choose  the  mapset  called  “PERMANENT.”  Every 
GRASS  location  has  a  mapset  called  PERMANENT,  which  is  reserved  to 
contain  any  final  maps.  You  will  be  asked  to  supply  information  such  as 
the  coordinate  system,  the  projection,  the  cell  resolution,  and  coordinates 
for  the  north,  south,  east,  and  west  extent  of  the  study  area.  You  can  find 
this  information  back  in  ArcCatalog  by  looking  at  the  spatial  attributes  of 
your  area  into  which  you  re-projected  your  maps  (e.g.,  DEM  or  NLCD). 

Enter  the  information  when  prompted:  Nad.83,  Data  Transformation  Pa¬ 
rameter:  6,  Zone _ ,  N,  then  enter  extents  a  little  big  for  the  area** 

It  is  very  important  to  set  the  resolution  to  30:  meters  and  ensure  that  the 
difference  between  the  N-S  and  E-W  need  to  be  equally  divisible  by  30. 

For  example,  if  the  N  is  36,000  and  the  S  is  30,000.  So  36,000  -  30,000  = 
6,000.  Also,  6,000  /  30  =  200.  Since  it  is  a  whole  number,  the  resolution 
is  30. 

Accept  the  region  when  prompted  and  then  exit  GRASS.  You  have  now 
successfully  set  up  a  new  GRASS  location! 

Import  the  ESRI  Data 

Start  GRASS  again  and  select  the  location  you  just  created  and  se¬ 
lect/create  a  new  mapset  that  you  must  name  “leamBase.”  By  convention, 
this  mapset  name  will  contain  all  of  the  base  maps  used  for  the  land  use 
change  analysis.  Also,  by  convention  all  map  names  will  be  specifically  as¬ 
signed  to  allow  for  various  scripts  to  work  properly.  Change  directory  to 
where  your  data  is  and  import  the  files  you  brought  in  using  r.in.shape  for 
Shapefiles  and  r.in.gdal  for  rasters  as  follows: 


v. in. shape  in=interstates . shp  out=inter states 

v. in. shape  in=otherroads . shp  out=otherroads  cat=CPCC3 

r.in.shape  in=nogrowth. shp  out=nogrowth 

r.in.shape  in=county.shp  out=county 

v. in. shape  in=boundary. shp  out=boundary 

r.in.shape  in=boundary. shp  out=boundary 

r.in.shape  in=mnbnd  out=mibnd 

r.in.gdal  in=landcover  out=landcover 

r.in.gdal  in=den  out=dem 
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These  files  will  be  stored  in  the  leamBase  mapset  folder  for  the  current  lo¬ 
cation  in  GRASS.  The  files  can  be  listed  using  “g.list  rast”  and  “g.list  vect”, 
assuming  “leamBase”  is  in  your  mapset  list  (see  “g.mapsets  help”).  Use  the 
GRASS  display  commands  (e.g.,  d.mon,  d.rast,  d.what.rast,  d.vect,  and 
d.what.vect)  to  check  your  results. 


Create  the  “City”  Files 

Required  analysis  input  maps  include  four  “city”  maps,  each  identifying 
attractor  areas  to  urban  development  based  on  “city  size.”  These  attractors 
are  points  of  locally  highest  residential  density  based  on  the  NLCD  land- 
cover  maps  and  discovered  with  the  following  GIS  analysis  script. 


#  This  creates  a  surface  of  "connectedness  to  local  residential" 

#  The  filter  looks  up  to  15  cells  away  and  assigns  an  importance  of  residential  depending 

#  on  distance  and  type  (e.g.,  the  1  or  2  value  from  above) .  It  multiplies  the  two  values 

#  for  each  cell  together  and  adds  the  results  for  all  neighboring  cells, 
r.mfilter  input=tl  output=t2  filter=neighborhoodFilter 

####  Next,  find  the  center  of  the  higher  attracting  areas 

#  First,  find  the  highest  values 

r.mfilter  input=t2  output=t3  filter=peakFilter 
r.mapcalc  t4='if (t2>(t3+400)  &&  t2>0) ' 

#  Then,  find  the  center  of  each  area. 

#  -  Make  clumps 
r. clump  t4  out=t5 

#  -  Find  centroids 

r.  volume  data=t3  clump=t5  site=centroids 

#  -  convert  to  cells 

s.  to. rast  centroids  out=t6  size=3 


#  Get  original  mfilter  output  numbers  for  each  centroid 
r.mapcalc  t7='if (t6, t2, 0) ' 


#  Divide  into  4  sets 

r.mapcalc  t8=round\  (4\*t7/'r. describe  -r  t7  |  grep  thru  |  sed  -e  's/.*  //'A) 


#  set  to  "small,  medium,  large,  and  x-large" 
r.mapcalc  smalljcity='if (t8<2, l,null () ) ' 
r.mapcalc  medium_city= ' if  (t8=2, l,null  () )  ' 
r.mapcalc  large_city='if  (t8=3,  l,null  () )  ' 
r.mapcalc  xl_city='if  (t8=4,  l,null  () )  ' 

The  two  filters  used  in  the  script  follow. 

The  first  is  "neighborhoodFilter, "  which  is  used  to  perform  a  distance-weighted  analysis  of  all 
residential  areas  around  each  cell  using  the  r.mfilter  program. 


TITLE  Find  dense  urban  areas 

MATRIX  29 

00000000000111111100000000 

00000000111222222211100000 

00000011222333333322211000 

00000112233444444433221100 

00001223344555555544332210 

00012234455666666655443221 

00112344566677777666544321 


0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
10  0 
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DIVISOR  1 
TYPE  P 


The  second  is  peakFilter,  which  finds  the  average  value  in  a  neighborhood 
surrounding  each  cell. 


TITLE 

MATRIX 

0  0  0 

Find 

29 

0  0  0  0 

average 

0  0  0  0 

in  neighborhood 

11111110 
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Process  the  Roads  Shape  File 

Roads  are  used  in  the  LEAM  processing  to  compute  travel  times  and  to 
provide  context  information  in  output  display  maps.  To  accommodate  the 
roads  file  must  be  imported  into  GRASS  as  two  vector  files,  interstates  and 
other  roads  (non-interstates).  These  are  then  converted  to  raster  files  for 
later  travel-time  processing  with  the  following  commands: 

#  List  the  categories  associated  with  the  roads  file 
v. in. shape  interstates . shp  -d 

v. in. shape  other roads . shp  -d 

#  Imports  the  roads  file  into  GRASS 

v. in. shape  interstates . shp  out=interstates  att  -o 

#  Creates  vector  support  files  for  interstates 
v. support  interstates 

v. in. shape  other roads . shp  out=otherroads  att=CPCC3  -o 

#  Creates  vector  support  files  for  otherroads 
v. support  otherroads 

#  Converts  both  to  raster 
v.to.rast  interstates  out=inter states 
v.to.rast  otherroads  out=otherroads 

Status 


At  this  point  the  following  maps  should  exist  as  part  of  the  leamBase 
mapset: 


Map  name 

Type 

Description 

small_city 

Raster 

Small  City 

medium_city 

Raster 

Medium  City 

large_city 

Raster 

Large  City 

xl_city 

Raster 

Extra  Large  City 

nogrowth 

Raster 

No  Growth 

dem 

Raster 

DEM 

mnbnd 

Raster 

Municipal  Boundaries 

county 

Raster 

Study  Area  /  Counties 

landcover 

Raster 

NLCD 

boundary 

Raster 

Boundary  of  some  area  of  interest 

otherroads 

Raster 

Non-interstates 

interstates 

Raster 

Interstate  highways 

boundary 

Vector 

Boundary  of  some  area  of  interest 

otherroads 

Vector 

Non-interstates 

interstates 

Vector 

Interstate  highways 

REPORT  DOCUMENTATION  PAGE  Form  Approved 

OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

16-09-2006  Final 


4.  TITLE  AND  SUBTITLE 

LEAMram™:  Land  use  Evolution  and  impact  Assessment  Model  Residential  Attractiveness 
Model 


3.  DATES  COVERED  (From  -  To) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

James  D.  Westervelt  and  Joseph  Rank 


5d.  PROJECT  NUMBER 

SERDP 


5e.  TASK  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  Research  and  Development  Center  (ERDC) 
Construction  Engineering  Research  Laboratory  (CERL) 

PO  Box  9005, 

Champaign,  1L  61826-9005 


5f.  WORK  UNIT  NUMBER 

CS-1257 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

ERDC/CERL  TR-06-28 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

SERDP  Program  Office  Commander,  Fort  Bragg 

901  North  Stuart  Street,  Suite  303  Fort  Bragg,  NC  28303 

Arlington,  VA  22203 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

SERDP 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


14.  ABSTRACT 

This  document  describes  a  GIS  approach  that  can  quickly  and  inexpensively  identify  areas  around  military  installations  that  are  attractive 
to  residential  development  -  based  on  proposed  regional  investments  and  policies.  Analysis  results  can  then  be  used  to  inform  regional 
planning  and  to  help  identify  long-term  impacts  of  residential  development  on  an  installation’s  future  opportunity  to  train  and  test.  Pro¬ 
posed  state  and  local  investment  in  regional  planning  policies  and  projects  can  alter  the  expected  patterns  of  growth  and  should  be  evalu¬ 
ated  with  respect  to  their  short  and  long-term  impacts  on  the  ability  to  conduct  training  and  testing  at  the  installation.  For  demonstration 
purposes,  the  approach  is  demonstrated  for  a  multi-county  area  around  Fort  Bragg,  NC. 
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